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Abstract. We present distance measurements to 71 high redshift type la supemovae discovered during the first year of the 
5-year Supernova Legacy Survey (SNLS). These events were detected and their multi-color light-curves measured using the 
MegaPrime/MegaCam instrument at the Canada-France-Hawaii Telescope (CFHT), by repeatedly imaging four one-square 
degree fields in four bands. Follow-up spectroscopy was performed at the VLT, Gemini and Keck telescopes to confirm the 
nature of the supemovae and to measure their redshift. With this data set, we have built a Hubble diagram extending to z = 1, 
with all distance measurements involving at least two bands. Systematic uncertainties are evaluated making use of the multi- 
band photometry obtained at CFHT. Cosmological fits to this first year SNLS Hubble diagram give the following results : 
£1m = 0.263 ± 0.042 (stat) ± 0.032 (sys) for a flat ACDM model; and w = -1.023 + 0.090 (stat) ± 0.054 (sys) for a flat 
cosmology with constant equation of state w when combined with the constraint from the recent Sloan Digital Sky Survey 
measurement of baryon acoustic oscillations. 
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the European Southern Observatory using the Very Large Telescope 
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2003B-Q-9, and GS-2003B-Q-8) obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under a cooperative agreement with the NSF on 
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1. Introduction 

The discovery of the acceleration of the Universe stands as a 
major breakthrough of observational cosmology. Surveys of 
cosmologically distant Type la supernovae (SNe Ia: lRiess et al.l 
Il998t IPerlmutter et al.lfl999l) indicated the presence of a new, 
unaccounted-for "dark energy" that opposes the self-attraction 
of matter and causes the expansion of the Universe to accel- 
erate. When combined with indirect measurements using cos- 
mic microwave background (CMB) anisotropics, cosmic shear 
and studies of galaxy clusters, a cosmological world model has 
emerged that describes the Universe as flat, with about 70% of 
its energy c ontained in the form of this cosmic dark energy (see 
for example ISeliaketal 120051) . 

Current projects aim at directly probing the nature of the 
dark energy via a determination of its equation of state pa- 
rameter - the pressure to energy-density ratio - w = px/px, 
which also defines the time dependence of the dark energy den- 
sity: px ~ a ~3(i+H)^ w here a is the scale factor. Recen t con- 
straints on w jKnop et al . 2003: Tonry et al. 2003: Bar ris et alJ 
l2004tlRiess etal.ll2004 are consistent with a very wide range 
of Dark Energy models. Among them, the historical cosmolog- 
ical constant {w - -1) is 10'^" to 10^*^ smaller than plausible 
vacuum energies predicted by fundamental particle theories. It 
also cannot explain why matter and dark energy have compara- 
ble densities today. "Dynamical A" models have been proposed 
(quintessence, k-essence) based on speculative field models, 
and some predict values of w above -0.8 - significantly differ- 
ent from -1. Measuring the average value of w with a precision 
better than 0.1 will permit a discrimination between the null 
hypothesis (pure cosmological constant, w - -1) and some 
dynamical dark energy models. 

Improving significantly over current SN constraints on the 
dark energy requires a ten-fold larger sample (i.e. o(lOOO) at 
0.2 < z < 1., where w is best measured), in order to signifi- 
cantly improve on statistical errors but also, most importantly, 
on systematic uncertainties. The traditional method of measur- 
ing distances to SNe la involves different types of observations 
at about 10 different epochs spread over nearly 3 months: dis- 
covery via image subtraction, spectroscopic identification, and 
photometric follow-up, usually on several telescopes. Many ob- 
jects are lost or poorly measured in this process due to the ef- 
fects of inclement weather during the follow-up observations, 
and the analysis often subject to largely unknown systematic 
uncertainties due to the use of various instruments and tele- 
scopes. 
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The Supernova Legacy Survey (SNLS)' was designed to 
improve significantly over the traditional strategy as follows: 
1) discovery and photometric follow-up are performed with a 
wide field imager used in "rolling search" mode, where a given 
field is observed every third to fourth night as long as it remains 
visible; 2) service observing is exploited for both spectroscopy 
and imaging, reducing the impact of bad weather. Using a sin- 
gle imaging instrument to observe the same fields reduces pho- 
tometric systematic uncertainties; service observing optimizes 
both the yield of spectroscopic observing time, and the light- 
curve sampling. 

In this paper we report the progress made, and the cosmo- 
logical results obtained, from analyzing the first year of the 
SNLS. We present the data collected, the precision achieved 
both from improved statistics and better control of system- 
atics, and the potential of the project to further reduce and 
control systematic uncertainties on cosmological parameters. 
Section|2ldescribes the imaging and spectroscopic surveys and 
their current status. Sections|3]and|3present the data reduction 
and photometric calibration. The light-curve fitting method, the 
SNe samples and the cosmological analysis are discussed in 
Section|5] A comparison of the nearby and distant samples used 
in the cosmological analysis is performed in Section|S]and the 
systematic uncertainties are discussed in Section^ 

2. The Supernova Legacy Survey 

The Supernova Legacy Survey is comprised of two compo- 
nents: an imaging survey to detect SNe and monitor their light- 
curves, and a spectroscopic program to confirm the nature of 
the candidates and measure their redshift. 

2.1. The imaging survey 

The imaging is taken as part of the deep component of 
the CFHT Legacy Survey dCFHTLSl l2002l) using t he one 
square degree imager, MegaCam JBoulad e et alJEo03h . In to- 
tal, CFHTLS has been allocated 474 nights over 5 years and 
consists of 3 surveys: a very wide shallow survey (1300 square 
degrees), a wide survey (120 square degrees) and a deep sur- 
vey (4 square degrees). The 4 pointings of the deep survey are 
evenly distributed in right ascension (Table The observa- 
tions for the deep survey are sequenced in a way suitable for 
detecting supernovae and measuring their light-curves: in every 
lunation in which a field is visible, it is imaged at five equally 
spaced epochs during a MegaCam run (which lasts about 18 
nights). Observations are taken in a combination of r^, im plus 
gM otzm filters (the MegaCam filter set; see Section|3 depend- 
ing on the phase of the moon. Each field is observed for 5 to 7 
consecutive lunations. Epochs lost to weather on any one night 
remain in the queue until the next clear observing opportunity, 
or until a new observation in the same filter is scheduled. 

During the first year of the survey, the observing efficiency 
was lower than expected and the nominal observation plan 
could not always be fulfilled. The scheduled im exposures (3 
X 3600 s plus 2 X 1800s per lunation) and exposures (5 

' see http://cfht.hawaii.edu/SNLS/ 
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Field 


RA(2000) 


Dec (2000) 


E(B-V) (MW) 


Dl 


02:26:00.00 


-04:30:00.0 


0.027 


D2 


10:00:28.60 


+02:12:21.0 


0.018 


D3 


14:19:28.01 


+52:40:41.0 


0.010 


D4 


22:15:31.67 


-17:44:05.0 


0.027 



Table 1 Coordinates and average Milky Way extinction (from 
ISchlegel et al.ll 19981) of fields observed by the Deep/SN com- 
ponent of the CFHTLS. 

epochs X 1500 s) were usually acquired. Assigned a lower pri- 
ority, gM and zm received less time than originally planned: on 
average only 2.2 epochs of 1050 s were collected per lunation 
in gM, and 2 epochs of 2700 s in zm; for the latter, the average 
ignores the D2 field and the D3 field in 2003, for which only 
fragmentary observations were obtained in zm- With efficiency 
ramping up, gM and zm approached their nominal rate in May 
2004, and since then the nominal observation plan (detailed in 
ISuUivan et al.l2005l) is usually completed. 

Observations and real-time pre-processing are performed 
by the CFHT s taff using the Elixir reduction pipeline 
jMagnier & Cuillandre 2004.) . with the data products immedi- 
ately available to the SN search teams. We have set up two inde- 
pendent real-time pipelines which analyze these pre-processed 
images. The detection of new candidates is performed by sub- 
tracting a "past" image to the current images, where the past- 
image is constructed by stacking previous observations of the 
same field. The key element of these pipelines is matching the 
point spread function of a new ex posure to the past-image. This 
is don e using the Alard algorithm jAlard & Lupto nPl998l lAiar3 
I2OO0I) for one of the pipelines, and using a non-parametric ap- 
proach for the other New candidates are detected and mea- 
sured on the subtraction images; detections are matched to 
other detections in the field, if any. One of the pipelines pro- 
cesses all bands on an equal footing, the other detects in the im 
band (which is deep enough for trigger purposes) and measures 
fluxes in the other bands. The two candidate lists are merged af- 
ter each epoch and typically have an overlap greater than 90% 
for iM(AB) < 24.0 after two epochs in a dark run. The reasons 
for one candidate being found by only one pipeline are usually 
traced to different masking strategies or different handling of 
the CCD overlap regions. 

2.2. Spectroscopic follow-up 

Spectroscopy is vital in order to obtain SN redshifts, and to 
determine the nature of each SN candidate. This requires ob- 
servations on 8-10 meter class telescopes due to the faintness 
of these distant supernovae. Spectroscopic follow-up time for 
the candidates presented in this paper was obtained at a va- 
riety of telescopes during the Spring and Fall semesters of 
2003 and the Spring semester of 2004. The principle spectro- 
scopic allocations were at the European Southern Observatory 
Very Large Telescope (program ID <171.A-0486>; 60 hours 
per semester), and at Gemini-North and South (Program-IDs: 
GN-2004A-Q-19, GS-2004A-Q-11, GN-2003B-Q-9, and GS- 
2003B-Q-8; 60 hours per semester). Spectroscopic time was 
also obtained at Keck-I and Keck-II (3 nights during each 



Spring semester) as the D3 field cannot be seen by VLT or 
Gemini-South. Further complementary spectroscopic follow- 
up observations were also obtained at Keck-I (4 nights in each 
of 2003 A, 2003B and 2004A) as part of a detailed s tudv of the 
intermediate redshift SNe in our sample jEUis et al.[ in prep.). 

Most of the observations are performed in long-slit mode. 
The detailed spectroscop ic classification of t hes e candidates is 
discussed elsewhere (see lHowell et al.lf2005l and lBasa et all in 
prep.). In summary , we consider two classes of events (see 
iHowell et al.l l200a for the exact definitions): secure SNe la 
events ("SN la"), and probable la events ("SN la*"), for which 
the spectrum matches a SN la better than any other type, but 
does not completely rule out other possible interpretations. All 
other events which were not spectroscopically identified as 
SN la or SN la* were ignored in this analysis. 

The imaging survey still delivers more variable candidates 
than can actually be observed spectroscopically. Hence, an ac- 
curate ranking of these candidates for further observations is 
essential. This ranking is performed to optimize the SN la yield 
of our allocations. Our method uses both a photometric se- 
lection tool (discussed in lSuUivan et aljE005l) which performs 
real-time light-curve fits to reduce the contamination of core- 
collapse SNe, and a database of every variable object ever de- 
tected by our pipelines to remove AGN and variable stars which 
are seen to vary repeatedly in long-timescale data sets (more 
than one year). 

SN la candidates fainter than im = 24.5 (likely at z > 
1) and those with very low percentage increases over their 
host galaxies (whe re identification is extremely difficult - see 
iHowell et alJ2005l) are usually not observed. With the real-time 
light-curve fit technique, approximately 70% of our candidates 
turned out to be SNe la. The possible biases asso ciated with 
this selection were studied in ISulUvan et all tOO^ and found 
to be negligible. 

2.3. The first year data set 

The imaging survey started in August 2003 after a few months 
of MegaCam commissioning. (Some SN candidates presented 
here were detected during the commissioning period.) This pa- 
per considers candidates with maximum light up to July 15* 
2004, corresponding approximatively to a full year of opera- 
tion. During this time frame, which includes the ramping-up 
period of the CFHTLS, about 400 transients were detected, 142 
spectra were acquired: 20 events were identified as Type II su- 
pernovae, 9 as AGN/QSO, 4 as SN Ib/c, and 91 events were 
classified as SN la or SN la*. The 18 remaining events have 
inconclusive spectra. Table0gives the 91 objects identified as 
SN la or SN la* during our first year of operation. 

3. Data reduction 

3.1. Image preprocessing 

At the end of each MegaCam run, the i mages are pre-processed 
agaiii at CFHT using the Elixir pipeline jMagnier & Cuillandrd 
I2OO4I) . This differs from the real-time reduction process de- 
scribed in Section 12.11 in that master flat-field images and 
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fringe-coiTection frames are constructed from all available data 
from the entire MegaCam run (including PI data). The Elixir 
process consists of flat-fielding and fringe subtraction, with an 
approximate astrometric solution also derived. Elixir provides 
reduced data which has a uniform photometric response across 
the mosaic (at the expense of a non-uniform sky background). 
This "photometric flat-field" correction is constructed using ex- 
posures with large dithers obtained on dense stellar fields. 

The SNLS pipelines then associate a weight map with each 
Elixir-processed image (i.e. each CCD from a given exposure) 
from the flat-field frames and the sky background variations. 
Bad pixels (as identifi ed by Elixir), cosmic rays (detected using 
the Laplacian filter of Ivan DokkumlEoOll) . satellite trails, and 
saturated areas are set to zero in the weigh t maps. An object 
catalo g is then produced using SExtractor jBertin & ArnoutsI 
11994 . and point-like objects are used to derive an image 
quality (IQ) estimate. The sky background map computed 
by SExtractor is then subtracted from the image. We addi- 
tionally perform aperture photometry on the objects of the 
SExtractor catalog for the purpose of photometric calibration 
(see Section]?}. 

3.2. Measurement of supernova fluxes 

For each supernova candidate, the image with the best IQ (sub- 
sequently called "reference") is identified, and all other images 
(both science images and their weight maps) are resampled to 
the pixel grid defined by this reference. The variations of the 
Jacobian of the geometrical transformations, which translate 
into photometric non-uniformities in the re-sampled images, 
are sufficiently small (below the millimag level) to be ignored. 
We then derive the convolution kernels that would match the 
PSF (modeled using the DAOPHOT package lstetsoiill 19871) of 
the reference image to the PSF of the other resampled science 
images, but we do not perform the convolutions. These con- 
volution kernels not only match the PSFs, but also contain the 
photometric ratios of each image to the reference. We ensure 
that these photometric ratios are spatially uniform by imposing 
a spatially uniform kernel integral, but allow for spatial ker- 
nel shape variati ons as t he im ages may have spatially varying 
PSFs. Foflowing lAlardI j2000l) . the kernel is fit on several hun- 
dred objects selected for their high, though unsaturated, peak 
flux. The kernel fit is made more robust by excluding objects 
with large residuals and iterating. 

Our approach to the differential flux measurement of a 
SN is to simultaneously fit all images in a given filter with 
a model that includes (i) a spatially variable galaxy (constant 
with time), and (ii) a time-variable point source (the super- 
nova). The model is described in detail in iFabbrol j200ll) . The 
shape of the galaxy and positions of both galaxy and supernova 
are fit globally. The intensity D, in a pixel p of image / is 
modeled as: 

Di,p^[{fiPref+g)®k]^+bi (1) 

where /) are the supernova fluxes, P„/ is the PSF of the refer- 
ence image centered on the SN position; ki is the convolution 
kernel that matches the PSF of the reference image to the PSF 



of image /; g is the intensity of the host galaxy in the refer- 
ence image, and bt is a local (sky) background in image The 
non parametric galaxy "model" g is made of independent pix- 
els which represent the galaxy in the best IQ image. All fluxes 
(/) are expressed in units of the reference image flux. 

The fit parameters are: the supernova position and the 
galaxy pixel values (common to all images), the supernova 
fluxes, and a constant sky background (different for each im- 
age). In some images in the series, the supernova flux is known 
to be absent or negligible; these frames enter the fit as "zero 
flux images" and are thus used to determine the values of the 
galaxy pixels. The least-squares photometric fit minimizes: 

^Yj^^'^p^^'-p-^'-p^^ (2) 

i.p 

where 1^, and W;,^ are the image and weight values of pixel p 
in image i, and the sums run over all images that contain the 
SN position, and all pixels in the fitted stamp of this image. 

Note that this method does not involve any real image con- 
volution: the fitted model possesses the PSF of the reference 
image, and it is the model that is convolved to match the PSF 
of every other image. We typically fit 50x50 galaxy pixels and 
several hundred images, and each SN fit usually has 2000 to 
3000 parameters. The fit is run once, 5cr outlier pixels are re- 
moved, and the fit is run again. 

The photometric fit yields values of the fit parameters along 
with a covariance matrix. There are obvious correlations be- 
tween SN fluxes and galaxy brightness, between these two pa- 
rameters and the background level, and between the SN posi- 
tion and the flux, for any given image. More importantly, the 
uncertainty in the SN position and the galaxy brightness intro- 
duces correlations between fluxes at different epochs that have 
to be taken into account when analyzing the light-curves. Note 
that flux variances and the coiTelations between fluxes decrease 
when adding more "zero flux images" into the fit. It will there- 
fore be possible to derive an improved photometry for most of 
the events presented in this paper, when the fields are observed 
again and more images without SN light are available. 

3.3. Flux uncertainties 

Once the photometric fit has converged, the parameter covari- 
ance matrix (including flux variances and covariances) is de- 
rived. This Section addresses the accuracy of these uncertain- 
ties, in particular the flux variances and covariances, which are 
used as inputs to the subsequent light-curve fit. 

The normalization of the parameter covariance matrix di- 
rectly reflects the normalization of image weights. We checked 
that the weights are on average properly normalized because 
the minimum per degree of freedom is very close to 1 (we 
find 1.05 on average). However, this does not imply mathe- 
matically that the flux uncertainties are properly normalized, 
because equation (|2} neglects the correlations between neigh- 
boring pixels introduced by image re-sampling. We considered 
accounting for these correlations; however, this would make 
the fitting code intolerably slow, as the resulting would be 
non-diagonal. Using approximate eiTors in least squares (such 
as ignoring correlations) increases the actual variance of the 
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estimators, but in the case considered here, the loss in photo- 
metric accuracy is below 1%. The real drawback of ignoring 
pixel correlations is that parameter uncertainties extracted from 
the fit are underestimated (since pixel correlations are posi- 
tive); this is a product of any photometry method that assumes 
uncorrelated pixels on re-sampled or convolved images. Our 
geometric alignment technique, used to align images prior to 
the flux measurement as described in Section IT2I uses a 3x3 
pixel quadratic re-sampling kernel, which produces output pix- 
els with an average variance of 80% of the input pixel variance, 
where the remaining 20% contributes to covariance in nearby 
pixels. We checked that flux variances (and covariances) com- 
puted assuming independent pixels are also underestimated by 
the same amount: on average, a 25% increase is required. 

In order to derive accurate uncertainties, we used the fact 
that for each epoch, several images are available which mea- 
sure the same object flux. Estimating fluxes on individual ex- 
posures rather than on stacks per night preserves the photomet- 
ric precision since a common position is fit using all images. 
It also allows a check on the consistency of fluxes measured 
within a night. We therefore fit a common flux per night to the 
fluxes measured on each individual image by minimizing a xl 
(where n stands for nights); this matrix is non-diagonal because 
the differential photometry produces correlated fluxes. The;;^^ 
contribution of every individual image is evaluated, and out- 
liers > 5cr (due to, for example, unidentified cosmic rays) are 
discarded; this cut eliminates 1 .4% of the measurements on av- 
erage. The covariance of the per-night fluxes is then extracted, 
and normalized so that the minimum per degree of freedom 
is 1. This translates into an "effective" flux uncertainty derived 
from the scatter of repeated observations rather than from first 
principles. If the only source of noise (beyond photon statistics) 
were pixel correlations introduced by image resampling, we 
would expect an average xl/Nd,,/ of 1.25, as all flux variances 
are on average under-estimated by 25%. Our average value is 
1.55; hence we conclude that our photometric uncertainties are 
only ~ 12% ( V(l-55/1.25) - 1) larger than photon statistics, 
leaving little margin for drastic improvement. 

Table |2 summarizes the statistics of the differential pho- 
tometry fits in each filter. The larger values of xl INdoj in im 
and zm probably indicate contributions from residual fringes. 
Examples of SNe la light-curves points are presented in 
Figures n and 121 showing SNe at z = 0.358 and z - 0.91 re- 
spectively. Also shown on these Figures are the results of the 
light-curves fits described in Section lSJl 

The next Section discusses how accurately the SN fluxes 
can be extracted from the science frames relative to nearby 
field stars, i.e. how well the method assigns magnitudes to SNe, 
given magnitudes of the field stars which are used for photo- 
metric calibration, called tertiary standards hereafter. 




3100 3150 3200 

JD 2450000+ 

Fig. 1 Observed light-curves points of the SN la SNLS-04D3fk 
in gM, I'm, Im and zm bands, along with the multi-color light- 
curve model (described in Section lSJl . Note the regular sam- 
pling of the observations both before and after maximum light. 
With a SN redshift of 0.358, the four measured pass-bands lie 
in the wavelength range of the light-curve model, defined by 
rest-frame U to R bands, and all light-curves points are there- 
fore fitted simultaneously with only four free parameters (pho- 
tometric normalization, date of maximum, a stretch and a color 
parameter). 




JD 2450000+ 

Fig. 2 Observed Ught-curves points of the SN la SNLS-04D3gx 
at z=0.91. With a SN redshift of 0.91, only two of the mea- 
sured pass-bands lie in the wavelength range of the light-curve 
model, defined by rest-frame U to R bands, and are therefore 
used in the fit (shown as solid lines). Note the excellent quality 
of the photometry at this high redshift value. Note also the clear 
signal observed in and even in gM, which correspond to cen- 
tral wavelength of respectively A ~ 3200A and A ~ 2500A in 
the SN rest-frame. 



3.4. Photometric alignment of supernovae relative to 
tertiary standards 

The SN flux measurement technique of Section ll!2l deIivers SN 
fluxes on the same photometric scale as the reference image. In 
this Section, we discuss how we measure ratios of the SN fluxes 



to those of the tertiary standards (namely stars in the SNLS 
fields). The absolute flux calibration of the tertiary standards 
themselves is discussed in Section|3 

The image model that we use to measure the SN fluxes (eq. 
^ can also be adapted to fit the tertiary standards by setting 
the "underlying galaxy model" to zero. We measure the fluxes 
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Band 


Average nb. 


Average nb. 


xl 


Central 




of images 


of epochs 


per d. o. f. 


wavelength 


8m 


40 


9.8 


1.50 


4860 




75 


14.4 


1.40 


6227 


Im 


100 


14.8 


1.63 


7618 


Zm 


60 


7.9 


1.70 


8823 



Table 2 Average number of images and nights per band for 
each SNLS light-curve. Note that there is less data in gM and 
Zm- The;t'^ column refers to the last fit that imposes equal fluxes 
on a given night. The expected value is 1.25 (due to pixel cor- 
relations), so we face a moderate scatter excess of about 12% 
over photon statistics. The larger values in Im and zm indicate 
that fringes play a role in this excess. The last column displays 
the average wavelength of the effective filters in A 

of field stars by running the same simultaneous fit to the im- 
ages used for the supernovae, but without the "zero-flux" im- 
ages, and without an underlying galaxy. As this fitting tech- 
nique matches that used for the SNe as closely as possible, 
most of the systematics involved (such as astrometric align- 
ment residuals, PSF model uncertainties, and the convolution 
kernel modeUng) cancel in the flux ratios. 

For each tertiary standard (around 50 per CCD), we ob- 
tain one flux for each image (as done for the SNe), expressed 
in the same units. From the magnitudes of these fitted stars, 
we can extract a photometric zero point for the PSF photom- 
etry for every star on every image, which should be identical 
within measurement uncertainties. Several systematic checks 
were performed to search for trends in the fitted zero-points as 
a function of several variables (including image number, star 
magnitude, and star color); no significant trends were detected. 
As zero-points are obtained from single measurements on sin- 
gle images, the individual measurements are both numerous 
and noisy, with a typical r.m.s of 0.03 mag; however since they 
have the same expectation value, we averaged them using a ro- 
bust fit to the distribution peak to obtain a single zero-point per 
observed filter 

To test how accurately the ratio of SN flux to tertiary stan- 
dard stars is retrieved by our technique, we tested the method 
on simulated SNe. For each artificial supernova, we selected a 
random host galaxy, a neighboring bright star (the model star), 
and a down-scale ratio (r). For half of the images that enter the 
fit, we superimposed a scaled-down copy (by a factor r) of the 
model on the host galaxy. We rounded the artificial position at 
an integer pixel offset from the model star to avoid re-sampling. 
We then performed the full SN fit (i.e. one that allows for an 
underlying galaxy model and "zero flux images") at the posi- 
tion of the artificial object, and performed the calibration star 
fit (i.e. one with no galaxy mode and no "zero-flux images") at 
the original position of the model star. This matches exactly the 
technique used for the measurement and calibration of a real 
SN. We then compared the recovered flux ratio to the (known) 
down-scale ratio. 

We found no significant bias as a function of SN flux or 
galaxy brightness at the level of 1 %, except at signal-to-noise 
(S/N) ratios (integrated over the whole light-curve) below 10. 
At a S/N ratio of 10, fluxes are on average underestimated by 



less than 1%; this bias rises to about 3% at a S/N ratio of 7. 
This small flux bias disappears when the fitted object position is 
fixed, as expected because the fit is then linear For this reason, 
when fitting zm light-curves of objects at z > 0.7, for which the 
S/N is expected to be low, we use the fixed SN position from 
that obtained from the Im and vm fits. 

Given the statistics of our simulations, the systematic un- 
certainty of SN fluxes due to the photometric method employed 
is less than 1 % across the range of S/N we encounter in real 
data, and the observed scatter of the retrieved "fake SNe" fluxes 
behaves in the same way as that for real SNe. Over a limited 
range of S/N (more than 100 integrated over the whole light- 
curve), we can exclude biases at the 0.002 mag level. Our upper 
limits for a flux bias have a negligible impact on the cosmologi- 
cal conclusions drawn from the sample described here, and will 
likely be improved with further detailed simulations. 

4. Photometric calibration 

The supernova light-curves produced by the techniques de- 
scribed in Section 13.21 are calibrated relative to nearby field 
stars (the tertiary standards). Our next step is to place these in- 
strumental fluxes onto a photometric calibration system using 
observations of stars of known magnitudes. 

4.1. Photometric calibration of tertiary standards 

Several standard star caU bration catalogs are av ailable in 
the literature, such as the iLandolJ (Il983[ Il992b|) Johnson- 
Cousin s (Vega-based) UBVRI system, or the ISmith et all 
il2002h u'g'r'i'z' AB-magnitude system which is used to cal- 
ibrate the Sloan Digital Sky Survey (SDSS). However, there 
are systema t ic erro rs affecting the transformations between the 
ISmifli et all ( l2002l) s ystem and the widely used Landolt sys- 
tem. As discussed in lFukugita et alJ lll996l) . these errors arise 
from various sources, for example uncertainties in the cross- 
calibration of the spectral energy distributions of the AB funda- 
mental standard stars relative to that of Vega. Since the nearby 
SNe used in our cosmological fits were extracted from the 
literature a nd are typically c alibrated using the standard star 
catalogs of ILandoltl ill 992bl) . we adopted the same calibra- 
tion source for our high-redshift sample. This avoids introduc- 
ing additional systematic uncertainties between the distant and 
nearby SN fluxes, which are used to determine the cosmologi- 
cal parameters. To eliminate uncertainties associated with color 
corrections, we derive magnitudes in the natural MegaCam fil- 
ter system. 

Both standard and science fields were repeatedly observed 
over a period of about 18 months. Photometric nights were 
selected using the CFHT "Skyprobe" instrument ( ICuillandr3 
l2003h . which monitors atmospheric transparency in the direc- 
tion that the telescope is pointing. Only the 50% of nights with 
the smallest scatter in transparency were considered. For each 
night, stars were selected in the science fields and their aper- 
ture fluxes measured and corrected to an airmass of 1 using the 
average atmospheric extinction of Mauna Kea. These aperture 
fluxes were then averaged, allowing for photometric ratios be- 
tween exposures. Stable observing conditions were indicated 
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Fig. 3 The calibration residuals — i.e. the residuals around the 
mean magnitude of each Deep field tertiary standard — in the 
bands gM, ru, iu and zm, for all CCDs and fields, with one 
entry per star and epoch. The dispersion is below 1% in gM, ru 
and iM, and about 1.5% in zm- 

by a very small scatter in these photometric ratios (typically 
0.2%); again the averaging was robust, with 5-cr deviations re- 
jected. Observations of the Landolt standard star fields were 
processed in the same manner, though their fluxes were not 
averaged. The apertures were chosen sufliciently large (about 
6" in diameter) to bring the variations of aperture corrections 
across the mosaic below 0.005 mag. However, since fluxes are 
measured in the same way and in the same apertures in science 
images and standard star fields, we did not apply any aperture 
correction. 

Using standard star observations, we first determined zero- 
points by fitting linear color transformations and zero-points to 
each night and filter, however with color slopes common to all 
nights. In order to account for possible non-linearities in the 
Landolt to MegaCam color relations, the observed color-color 
relations were then compared to synthetic ones derived from 
spectrophotometric standards. This led to shifts of roughly 0.01 
in all bands other than gM, for which the shift was 0.03 due to 
the nontrivial relation to B and V. 

We then applied the zero-points appropriate for each night 
to the catalog of science field stars of that same night. These 
magnitudes were averaged robustly, rejecting 5-cr outliers, and 
the average standard star observations were merged. Figure |3] 
shows the dispersion of the calibration residuals in the gM, ru, 
iM and zm bands. The observed standard deviation, which sets 
the upper bound to the repeatability of the photometric mea- 
surements, is about or below 0.01 mag in gM, fM and im, and 
about 0.016 mag in zm- 

For each of the four SNLS fields, a catalog of tertiary 
standards was produced using the procedure described above. 



These catalogs were then used to calibrate the supernova fluxes, 
as described in Section 13.41 The dominant uncertainty in the 
photometric scale of these catalogs comes from the determina- 
tion of the color-color relations of the standard star measure- 
ments. For the gM, rM and Im bands, a zero-point offset of 0.01 
mag would easily be detected; hence we took this value as a 
conservative uncertainty estimate. The zm band is aff'ected by a 
larger measurement noise, and it is calibrated with respect to / 
and R - I Landolt measurements. We therefore attributed to it 
a larger zero point uncertainty of 0.03 mag. 

The MegaCam shutter is designed to preserve the mosaic 
illumination uniformity. Nevertheless, the shutter precision is a 
potential source of systematic uncertainties, given (1) the pos- 
sible non uniformities due to the shutter motion and (2) the ex- 
posure time differences between the calibration images (a few 
seconds) and the science images (hundreds of seconds). For 
MegaCam, the actual exposure time is measured and reported 
for each exposure, using dedicated sensors . The shutter preci- 
sion was investigated by ICuillandr3 ( l2005l) and it was shown 
that the non-uniformity due to the shutter is less than 0.3% 
across the mosaic. Short and long exposures of the same fields 
were also compared. The systematic flux diff'erences between 
the exposures were found to be below 1% (r.m.s). 

4.2. The MegaCam and Landolt instrumental filters 

As the supernova fluxes are measured in the instrumental fil- 
ter system, the MegaCam transmission functions (up to an ar- 
bitrary constant) are needed in order to correctly interpret the 
SN photometry. Similarly, for the published nearby supernovae 
which are reported in Landolt magnitudes, the filter responses 
of the Landolt system are required. 

For the MegaCam filters, we used the measurements pro- 
vided by the manufacturer, multiplied by the CCD quantum 
efficiency, the MegaPrime wide-field corrector transmission 
function, the CFHT primary mirror reflectivity, and the aver- 
age atmospheric transmission at Mauna Kea. As an additional 
check, we computed synthetic MegaCam-SDSS color terms us- 
ing the synthe tic transmissions of the SDSS 2.5-m telescope 
llSDSSll2004hl) and spectrophotometric stan dards taken from 
IPicklesI ( 11993) and 'Oun n & StrvkeJ (Il983l). Since the SDSS 
science catalog ( Finkbei neTet al.ll2004 iRaddick 2001 ISDSSi 
l2004j) shares thousands of objects with two of the four fields 
repeatedly observed with MegaCam, we were able to compare 
these synthetic color transformations with the observed trans- 
formations. We found a good agreement, with uncertainties at 
the 1% level. This constrains the central wavelengths of the 
MegaCam band passes to within 10 to 15 A with respect to the 
SDSS 2.5m band passes. 

The choice of filter band passes to use for Landolt-based 
observations is not unique. Most previous sup ernova cosmol- 
ogy works assumed that the determinations of lBesselH (Il990l) 
describe the effective Landolt system well, although the author 
himself questions this fact, explicitly warning that the Landolt 
system "is not a good match to the standard syst em" — i.e. the 
historical Johnsons-Cousins system. Fortunatelv.^Ha muv et al.l 
(1I992.. .19940 provide spectrophotometric measurements of a 
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few objects measured in iLandold (Il992ah : this enabled us to 
compare synthetic magnitudes computed using Bessell trans- 
missions with Landolt measurements of the same objects. This 
comparison reveals small residual color terms which vanish if 
the B, V, R and / Bessell filters are blue-shifted by 41, 27, 
21 and 25 A respectively. Furthermore, if one were to assume 
that the Bessell filters describe the Landolt system, this would 
lead to synthetic MegaCam-Landolt color terms significantly 
different from the measured ones; the blue shifts determined 
above bring them into excellent agreement. We therefore as- 
sumed that the Landolt catalog magnitudes refer to blue-shifted 
Bessell filters, with a typical central wavelength uncertainty of 
10 to 15 A, corresponding roughly to a 0.01 accuracy for the 
color terms. 

4.3. Converting magnitudes to fluxes 

Given the variations with time of the cosmological scale factor 
a{t), one can predict the evolution with redshift of the observed 
flux of classes of objects of reproducible luminosity though not 
necessarily known. This is why the cosmological conclusions 
that can be drawn from flux measurements rely on flux ratios 
of distant to nearby SNe, preferably measured in similar rest- 
frame pass-bands. The measured SNe magnitudes must there- 
fore be converted to fluxes at some point in the analysis. 

The flux in an imaginary rest-frame band of transmission 
Trest for a SN at redshift z is deduced from the magnitude 
m(Tgi,s) measured in an observer band of transmission T„hs via: 

j<PsN('l)Tres,U)dA f , ,^ (3) 



7 (l>sNWTohs{Ail+z))dA 



J <Pref{A)Tohs' 



{X)dA 



where (fisN is the spectrum of the SN, nirefiT) is the magnitude 
of some reference star that was used as a calibrator, and (p,-ef is 
its spectrum. In this expression, the product of the first and third 
terms gives the integrated flux in the observed band, and the 
second term scales this integrated flux to the rest-frame band. 
We measure only m{Tohs) — mrefiTots) (if the reference star is 
directly observed), or only m{Tohs) (if a non-observed star - 
e.g. Vega - is used as the reference). The reference spectrum, 
(pref, must be taken from the literature, as well as nirefiTobs) if 
the reference is not directly observed. The supernova spectrum, 
(psN, is taken to be a template spectrum appropriately warped 
to reproduce th e observed color of the SN (as described in 
IGuv et ai]l2005l) . The quantity /(z, Tresi) scales as the inverse 
square of a luminosity distance: 



f{Z\,Trest) IdiiZl) 



f{Z2,Trest) \'^l(Zi) 



(4) 



This conversion of a measured magnitude to a rest-frame 
flux (or a rest-frame magnitude) is usually integrated in the so- 
called cross-filter k-corrections /Ki m et al.ll996 ': 'Nugen t et"al] 
I2OO2I). In our c ase, it is integrated in the light-curve fit 
( IGuv et alJ2005l) . (See lGuv et aiT ( l2005l) for a discussion of the 
precise definitions of spectra and transmissions that enter into 

f{Z, Trest)-) 



Inspecting eq. |3 we first note that the normalizations of 
Tabs and (ps n cancel. The width of Tabs is a second order effect. 
When forming the ratio of two such quantities for two different 
SN, the normalization of (prej does not matter, nor the normal- 
ization of Trest, provided the same Trest is chosen for both ob- 
jects. The width of Trest matters only at the second order. The 
factors that do enter as first order effects are: 

- / <PreM)Tobs,l(A(l + Zl))dA/J (prefU)Tabs,2W + ZlWA , 

which requires both the spectrum of a reference and the 
band passes of the observing systems, i.e. to first order, their 
central wavelengths, 

- nirejiTohs.i) - nhefiTobs,!) , i-C- the color of the reference. 
When comparing distant and nearby SNe, we typically rely 
on B - R or B - I colors, 

- and obviously, the SNe measured magnitudes, or, more pre- 
cisely, their difference. 

We choose to use Vega as the refe rence star An accu- 
rate spectrum of Vega was assembled bv lHavesI ( Il985h . Some 
subtle differences are found by a more recent HST measure- 
ment ('Bohlin & Gillilandll2004 but they only marginally af- 
fect broadband photom e try: di fferences within the 1 % uncer- 
tainty quoted in iHavej ill 9851) are found and we will assign 
this uncertainty to the Vega broadband fluxes. We use the HST- 
based measurement because it extends into the UV and NIR 
and hence is safe for the blue side of the U band and in the zm 
band. For Vega, we adopt the magnitudes (U,B,y jirJr) - (0.02, 
0.03, 0.03, 0.03, 0.024) jFukugita et alJ ill 994 and references 
therein). For other bands, a simple interpolation is adequate. 
Note that only Vega colors impact on cosmological measure- 
ments. 

A possibl e shortcut consists in relyin g on spectrophotomet- 
ric standards llHamuv et al.ll 19921 1 1994 which also have mag- 
nitudes on the Landolt system llLandolt Il992ah . When we com- 
pare synthetic Vega magnitudes of these objects with the pho- 
tometric measurements, we find excellent matching of colors 
(at better than the 1 % level), indicating that choosing Vega or 
spectrophotometric fluxes as the reference makes little practi- 
cal difference. 



4.4. Photometric calibration summary 

We constructed catalogs of tertiary standard stars in the SNLS 
fields, expressed in MegaCam natural magnitudes, and defined 
on the Landolt standard system. The repeatability of measure- 
ments of a single star on a given epoch (including measurement 
noise) is about or below 0.01 mag rm.s in gM, fM and Im, and 
about 0.016 mag in zm- From standard star observations, we 
set conservative uncertainties of the overall scales of 0.01 mag 
in gM, fM and and 0.03 in z^. The MegaCam central wave- 
lengths are constrained by color terms with respect to both the 
SDSS 2.5m telescope and the Landolt catalog to within 10 to 
15 A. The central wavelengths of the band pass es of the Landol t 
catalog are found slightly offset with respect to lBessel] l ll99(]l) . 
using spectrophotometric measurements of a subsample of this 
catalog. 
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5. Light-curve fit and cosmological analysis 

To derive the brightness, light-curve shape and SN color es- 
timates required for the cosmological analysis, the time se- 
quence of photometric measurements for each SN was fit us- 
ing a SN light-curve model. This procedure is discussed in this 
section together with the nearby and distant SN la samples se- 
lection and the cosmological analysis. 

5.1. The SN la light-curve model 

We fit the SN la light-cur ves in two or mo re bands using the 
SALT light-curve model llGuv et aljEool which returns the 
supernova rest-frame B-band magnitude m*^, a single shape pa- 
rameter s and a single color parameter c. The supernova rest- 
frame B-band magnitude at the date of its maximum luminosity 
in B is defined as: 

- , , ( f(z, T*g, t = t,„ax,B 

irig = -2.5 logio r 

where T*g{A) = Tb(A/(1+z)) = Tre,tiB) is the rest-frame B-band 
transmission, and f{z, T^, t - t„,a.x,B) is defin ed by eg. El The 
stretch factor s is similar to that described in IPerlmutter et all 
lll997l) : it parameterizes the brighter-slower relation, originally 
described in Iphillips, (.IQQD . by stretching the time axis of a 
unique light-curve template; i = 1 is define d in rest-frame B 
for the time interval -15 to +35 days using the lGoldhaber et al] 
( l200lh B-band template. For bands other than B, stretch is a 
parameter that indexes light-curve shape variability. The rest- 
frame color c is defined by c = (Z? - V)Bma.K + 0.057: it is a 
color excess (or deficit) with respect to a fiducial SN la (for 
which B - V = -0.057 at B-band maximum light). Note that 
the color c is not just a measure of host galaxy extinction: c can 
accommodate both reddening by dust and any intrinsic color ef- 
fect dependent or not on s. The reference value (-0.057) can be 
changed without changing the cosmological conclusions, given 
the distance estimator we use (see Section|^3- 

The light-curve model was trained on very nearby super- 
novae (mostly a t z < 0.015) published in the literature (see 
iGuv et aljEool for the selection of these objects). Note that 
these training objects were not used in the Hubble diagram de- 
scribed in this paper. The SALT light-curve model generates 
light-curves in the observed bands at a given redshift, SALT 
also incorporates corr ections for the Milky W ay extinction, us- 
ing the dust map s of IS chlegel et al. (1999) coupled with the 
extinction law of ICardelli et al. Q2££)- The rest-frame cover- 
age of SALT extends from 3460 to 6500 A (i.e. slightly blue- 
wards from U to R). We require that photometry is available 
in at least 2 measured bands with central wavelengths within 
this wavelength range to consider a SN for the cosmologi- 
cal analysis. Light curves in the zm band become essential 
for z > 0.80, since at these redshifts, corresponds to rest- 
frame A < 3460 A. All observed bands are fitted simultane- 
ously, with common stretch and color parameters, global in- 
tensity and date of B-band maximum light. Making use of U-, 
B- and V-band mea surements of nearby SNe la from the lit- 
erature (mostly from'Hamuv et al."l99(j':' Riess et al.l 'l999^.'jha' 
.Guv et al.. L2005.) have constructed a distance estimator 



using either U- and B-band data or B- and V-band which shows 
a dispersion of 0.16 mag around the Hubble line. The fitted 
global intensity is then translated into a rest-frame-B observed 
magnitude at maximum light (nig) which does not include any 
correction for brighter-slower or brighter-bluer relations. 

The light-curve fit is carried out in two steps. The first fit 
uses all photometric data points to obtain a date of maximum 
light in the B-band. All points outside the range [-15, +35] 
rest-frame days from maximum are then rejected, and the data 
refit. This restriction avoids the dangers of comparing light- 
curve parameters derived from data with different phase cover- 
age: nearby SNe usually have photometric data after maximum 
light, but not always before maximum when the SN is rising, 
and almost never before -15 days. By contrast, SNLS objects 
have photometric sampling that is essentially independent of 
the phase of the light-curve because of the rolling-search ob- 
serving mode, though late-time data (in the exponential tail) 
often has a poor S/N, or is absent due to field visibility. 

5.2. The SN la samples 

The cosmological analysis requires assembling a sample of 
nearby and distant SNe la. 

We assembled a nearby SN la sample from the literature. 
Events with redshifts below z = 0.015 were rejected to limit 
the influence of peculiar velocities. We further retained only 
objects whose first photometric point was no more than 5 days 
after maximum light. To check for possible biases that this lat- 
ter procedure might have introduced, we fitted subsets of data 
from objects with pre-maximum photometry. Our distance es- 
timator (see Section 15. 4> was found to be unaffected if the 
data started up to 7 days after maximum light. A sample of 
44 nearby SNe la matched our requirements. Table|8]gives the 
SN name, redshift and filters used in the light-curve fits, as well 
as fitted rest-frame B-band magnitude and values of the param- 
eters s and c. 

For this paper, we considered only distant SNe la that 
were discovered and followed during the first year of SNLS 
since this data set already constitutes the largest well con- 
trolled homogeneous sample of distant SN la. As discussed 
in Section |231 91 objects were spectroscopically identified as 
"la" or "la*", with a date of maximum light before July 15, 
2004. Ten of these are not yet analyzed: 5 because images un- 
contaminated by SN light were not available at the time of this 
analysis, and 5 due to a limitation of our reduction pipeline 
which does not yet handle field regions observed with different 
CCDs. Six SNe have incomplete data due to either instrument 
failures, or persistent bad weather and two SNe, SNLS-03D3bb 
and SN LS-03D4ci. which happen to be spectroscopically pecu- 
liar (see lEUis et al.l in prep.) have photometric data incompati- 
ble with the light-curve model. 

The resulting fit parameters of the remaining 73 "Ia"+"Ia*" 
SNe are given in Table |9] and examples of light-curves mea- 
sured in the four MegaCam bands are shown in Figures ^ and 
121 together with the result of the light-curve fit. 
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5.3. Host galaxy extinction 

There is no consensus on how to correct for host galaxy ex- 
tinction affectin g high redshift SNe la. The pioneering SN cos- 
mology papers jRiess et aDll998i IPerlmutter et alJ ll999^ typi- 
cally observed in only one or two filters, and so had little or no 
color information with which to perform extinction corrections. 
Subsequent papers either selecte d low-extinction sub samples 
based on host galaxy diagnostics dSullivan et al.l20()3l) . or used 
multicolor information together with an assumed color of an 
unredde ned SN to rnake extinction correctio ns on a subset of 
the data dKnop et al.l2003UTonrv et alj2003h . 

These techniques have their drawbacks: the intrinsic color 
of SNe la has some dispersion, and measured colors often 
have large statistical errors in high-redshift data sets. When 
these two color uncertainties are multiplied by the ratio of to- 
tal to selective absorption, Rb - 4-, the resulting error can be 
very lar ge. To circumvent this, some studies i ised Bayesian pri- 
ors (e.g.kiess et alJll998t iTonrv et al.ll2003t iRiess et al.ll2004t 
iBarris et alJl2004. Other authors argue that this b iases the re- 
sults re. g. IPerlmutter et al.ll999HKnop et al.l2003l) . 

Here we employ a technique that makes use of color in- 
formation to empirically improve distance estimates to SNe la. 
We exploit the fact that the SN color acts in the same direc- 
tion as reddening due to dust - i.e. redder SNe are intrinsically 
dimmer, brighter SNe are intrinsically bluer jTripp & BranchI 
Il999l) . By treating the correction between color and brightness 
empirically, we avoid model-dependent assumptions that can 
both artificially inflate the errors and potentially lead to biases 
in the determination of cosmological parameters. Because we 
have more than one well-measured color for several SNe, we 
can perform consistency checks on this technique - distances 
from multiple colors should, and do, agree to a remarkable de- 
gree of precision (Section lOl . 



5.4. Cosmological fits 

From the fits to the light-curves (Section [5.1> . we computed a 
rest-frame-B magnitude, which, for perfect standard candles, 
should vary with redshift according to the luminosity distance. 
This rest-frame-B magnitude refers to ofoiervet/ brightness, and 
therefore does not ac count for bright er-slower and brighter- 
bluer correlations (see lGuv et al .120051 and references therein). 
As a distance estimator, we use: 

A'B — - M + a(s — 1) -ySc 

where m* , s and c are derived from the fit to the hght curves, 
and a, /? and the absolute magnitude M are parameters which 
are fitted by minimizing the residuals in the Hubble diagram. 
The cosmological fit is actually performed by minimizing: 



ob jects 



{piB-5\og,a{dL(B,z)llOpc)f 



where stands for the cosmological parameters that define the 
fitted model (with the exception of //o), di is the luminosity 
distance, and cr,„, is the intrinsic dispersion of SN absolute 
magnitudes. We minimize with respect to 9, a, /3 and M. Since 



di scales as I/Hq, only M depends on Ho- The definition of 
cr^iMB), the measurement variance, requires some care. First, 
one has to account for the full covariance matrix of nig, s and 
c from the light-curve fit. Second, o-(fj.B) depends on a and /3; 
minimizing with respect to them introduces a bias towards in- 
creasing errors in order to decrease the;;^^, as originally noted 
in lTripp| ( |l998l) . When minimizing, we therefore fix the values 
of a and yS entering the uncertainty calculation and update them 
iteratively. crifiB) also includes a peculiar velocity contribution 
of 300 km/s. cr,„, is introduced to account for the "intrinsic dis- 
persion" of SNe la. We perform a first fit with an initial value 
(typically 0.15 mag.), and then calculate the cr,„, required to 
obtain a reduced - 1. We then refit with this more accurate 
value. We fit 3 cosmologies to the data: a A cosmology (the pa- 
rameters being Qm and Qa), a flat A cosmology (with a single 
parameter Om), and a flat w cosmology, where w is the constant 
equation of state of dark energy (the parameters are Qm and w). 

The Hubble diagram of SNLS SNe and nearby data is 
shown in Figure |4] together with the best fit A cosmology for 
a flat Universe. Two events lie more than 3 cr away from the 
Hubble diagram fit: SNLS-03D4au is 0.5 mag fainter than the 
best-fit and SNLS-03D4bc is 0.8 mag fainter. Although, keep- 
ing or removing these SNe from the fit has a minor effect on 
the final result, they were not kept in the final cosmology fits 
(since they obviously depart from the rest of the population) 
which therefore make use of 44 nearby objects and 71 SNLS 
objects. 

The best-fitting values of a and /3 are a - 1.52 ± 0.14 
and y6 = 1.57 + 0.15, comparable with prev ious works u sing 




199Rft . As 
and refer- 



similar distance estimators (see f or example 
discussed by several authors (see IGuv et al\ 
ences therein), the value of /3 does differ considerably from 
Rb = 4-, the value expected if color were only affected by 
dust reddening. This discrepancy may be ar i indicator of intriii - 
sic color variations in the SN sample (e.g. iNobili et alJl2003l) . 
and/or variations in Rb- For the absolute magnitude M, we ob- 
tain M = -19.31 +0.03 + 5 logio/ivo. 

The parameters a, p and M are nuisance parameters in the 
cosmological fit, and their uncertainties must be accounted for 
in the cosmological error analysis. The resulting confidence 
contours are shown in Figures|5]and|6l together with the prod- 
uct of these confidence estimates with the probability distribu- 
tion from bary on acoustic oscillations (BAO) measured in the 
SDSS (Eq. 4 in lEisenstein et al.l2005l) . We impose w = -1 for 
the {Qm, ^a) contours, and = for the {D.m, w) contours. 
Note that the constraints from BAO and SNe la are quite com- 
plementary. The best-fitting cosmologies are given in Table|3 



fit 


parameters (stat only) 




(0.31 ±0.21,0.80 + 0.31) 




(-0.49 + 0.12,1.11 +0.52) 


i^lu, flat 


Qm = 0.263 + 0.037 


(f^M, ^a) + BAO 


(0.271 +0.020,0.751 +0.082) 


(nM,w)+BAO 


(0.271 +0.021,-1.023 + 0.087) 



Table 3 Cosmological parameters and statistical errors of 
Hubble diagram fits, with the BAO prior where applicable. 
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Fig. 4 Hubble diagram of SNLS and nearby SNe la, with var- 
ious cosmologies superimposed. The bottom plot shows the 
residuals for the best fit to a flat A cosmology. 




0.5 1 

Fig. 5 Contours at 68.3%, 95.5% and 99.7% confidence levels 
for the fit to an (Qm, ^^a) cosmology from the SNLS Hubble di- 
agram (solid contours) , the SDSS baryon acoustic oscillations 
( lEisenstein et aljE005l dotted lines), and the joint confidence 
contours (dashed lines). 



Using Monte Carlo realizations of our SN sample, we 
checked that our estimators of the cosmological parameters are 
unbiased (at the level of 0. 1 cr), and that the quoted uncertain- 
ties match the observed scatter. We also checked the field-to- 
field variation of the cosmological analysis. The four Qm val- 
ues (one for each field, assuming Qk = 0) are compatible at 
37% confidence level. We also fitted separately the la and la* 
SNLS samples and found results compatible at the 75% confi- 
dence level. 

We derive an intrinsic dispersion, cr„„ = 0.13 ± 0.02, ap 



preciably smaller than previously measure d (|Riess_eLalJ 
Perlmutter et"al]ll999t iTonrv et all l2003t iBarris et all 



1998 



2004 



Riess'etaLn2004 ). The intrinsic dispersions of nearby only 
(0.15±0.02) and SNLS only (0.12+0.02) events are statistically 
consistent although SNLS events show a bit less dispersion. 

A notable feature of Figurel^is that the error bars increase 
significantly beyond z=0.8, where the zm photometry is needed 



to measure rest-frame B - V colors. The zm data is afifected by 
a low signal-to-noise ratio because of low quantum efficiency 
and high sky background. For z > 0.8, cr{{B - V)resiframe) - 
1 .6cr{iM-ZM), because the lever arm between the central wave- 
lengths of iM and Zm is about 1 .6 times lower than for B and V. 
Furthermore, errors in rest-frame color are scaled by a further 
factor of p - 1.6 in the distance modulus estimate. With a typ- 
ical measurement uncertainty cr{zM) - 0.1, we have a distance 
modulus uncertainty cr(yu) > 0.25. Since the fall 2004 semester, 
we now acquire about three times more zm data than for the 
data in the current paper, and this will improve the accuracy of 
future cosmological analyses. 

The distance model we use is linear in stretch and color. 
Excluding events at z > 0.8, where the color uncertainty is 
larger than the natural color dispersion, we checked that adding 
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Fig. 6 Contours at 68.3%, 95.5% and 99.7% confidence lev- 
els for the fit to a flat (Qm,w) cosmology, from the SNLS 
Hubble diagram alone, from the SDSS baryon acoustic oscil- 
lations alone CEisenstein et al...2005J . and the joint confidence 
contours. 



quadratic terms in stretch or color to the distance estimator de- 
creases the minimum by less than 1 . We hence conclude that 
the linear distance estimator accurately describes our sample. 

Since the distance estimator we use depends on the color 
parameter c, residuals to the Hubble Diagram are statistically 
correlated to c. The correlation becomes very apparent when 
the c measurement uncertainty dominates the distance uncer- 
tainty budget, as happens in our sample when z > 0.8. We 
checked that the measurement uncertainties can account for the 
observed residual-c correlation at z > 0.8. Because of this cor- 
relation, color selected sub-samples mechanically lead to bi- 
ased estimations of cosmological parameters. 

6. Comparison of nearby and distant SN properties 

6.1. Stretch and color distributions 

The distributions of the shape and color parameters - s and c 
as defined in Section lsHI - are compared in Figures0and|8]for 
nearby objects and for SNLS supernovae at z < 0.8 for which 
c is accurately measured. These distributions look very similar, 
both in central value and shape. The average values for the two 
samples differ by about Icr in stretch and 1 .5cr in color: we find 
that distant supernovae are on average slightly bluer and slower 
than nearby ones. The statistical significance of the differences 
is low and the differences can easily be interpreted in terms of 
selection effects rather than evolution. The evolution of average 
s and c parameters with redshift is shown in Section fT^I stretch 
is not monotonic, and color seems to drift towards the blue with 
increasing redshift. We show in Section 17.41 that the bulk of 




stretch 

Fig. 7 The stretch (s parameter) distributions of nearby (hashed 
blue) and distant (thick black with filled symbols) SNLS SNe 
with z < 0.8. These distributions are very similar with averages 
of 0.920 + 0.018 and 0.945 ± 0.013, respectively (Icr apart). 




-0.4 -0.2 0.2 0.4 

color 



Fig. 8 The color (c parameter) distributions for nearby (hashed 
blue) and distant (thick black with filled symbols) SNe with 
z < 0.8. These distributions are very similar, with averages of 
0.059 + 0.014 and 0.029 + 0.015, respectively (1.5cr apart). 



this effect can be reproduced by selection effects applied to an 
unevolving population. 

6.2. Brighter-slower and brighter-bluer relationships 

Figures |5] and [TO] compare the nearby and distant samples in 
the stretch-magnitude and color-magnitude planes. There is no 
significant difference between these samples. 

In Figure m two of the SNLS events (SNLS-04Dlag and 
SNLS-04D3oe) have a color value, c, smaller than -0.1. These 
supernovae are both classified as secure la. There are no SNe la 
in the nearby sample that are this blue. Figure ^| shows that 
these events lie on the derived brighter-bluer relation. Although 
they are brighter than average, fitting with or without these two 
events changes the cosmological results by less than 0.1 cr. 

6.3. Compatibility of SN colors 

The measurement of distances to high redshift SNLS SNe in- 
volves the rest-frame U band. The MegaCam band shifts 
from rest-frame B at z=0.5 to rest-frame U at z=0.8. Within 
this redshift range, distances are estimated mainly using im and 
rM, the weight of zm being affected by high photometric noise; 
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Stretch 

Fig. 9 Residuals in the Hubble diagram as a function of stretch 
(s parameter), for nearby (blue open symbols) and distant 
(z < 0.8, black filled symbols). This diagram computes dis- 
tance modulus jUb without the stretch term a{s - 1), and returns 
the well-known brighter-slower relationship with a consistent 
behavior for nearby and distant SNe la. 




m 



-0.2 0.2 0.4 

color 

Fig. 10 Residuals in the Hubble diagram as a function of color 
(c parameter), for nearby (blue open symbols) and distant (z < 
0.8, black filled symbols). This diagram computes distance 
modulus jUb without the color term /3c, and returns the brighter- 
bluer relationship with a consistent behavior for nearby and dis- 
tant SNe la. Notice that the bluest SNLS objects are compatible 
with the bulk behavior. 



2. We fit the two reddest bands (BV for nearby objects), with 
the stretch, and date of maximum being fixed at the previ- 
ously obtained values. From the fitted light-curve model we 
extract the expected rest-frame U band magnitude at maxi- 
mum B light, Ubv- 

3. We fit the two bluest bands, {UB for nearby objects), still 
with the stretch and date of maximum fixed. From this fit, 
we extract the expected rest-frame U band magnitude at 
maximum B light. Since it matches the measurement when 
the actual U flux is measured, we call it U,„eas- 

The test quantity is AU-i = Ubv — U,neas, i-C- the "predicted" 
U (derived from B and V) minus the measured U brightness. 
Forcing both quantities to be measured with the same stretch 
and B maximum date is not essential, but narrows the distribu- 
tion of residuals. A residual of zero means that the three mea- 
sured bands agree with the light-curve model for a certain pa- 
rameter set, and hence that the distance estimate will be identi- 
cal for the two different color fits. 

There are 10 SNLS "intermediate" redshift events at 0.25 < 
z < 0.4, where gMrMiM sample the UBV rest-frame region, and 
17 "distant" events at 0.55 < z < 0.8, where UBV shifts to the 
rMiMiM triplet. We also have at our disposal a sample of 28 
"nearby" objects measured mUBV, both from the nearby sam- 
ple described in Table |8l and also from the light-curve model 
trai ning sample whic h consists mainly of very nearby objects 
(see lGuv et al.*200"5'). Figure^Jdisplays the value of AC/3 as a 
function of redshift and Table |4] summarizes the averages and 
dispersions. A very small scatter (about 0.033) is found for the 
intermediate redshift sample. The nearby and distant samples 
exhibit larger scatters; the nearby sample is probably affected 
by the practical difficulties in calibrating U observations, and 
our distant sample is affected by the poor S/N in the zm band. 
We conclude from this study that our light curves model ac- 
curately describes the relations between the supernovae colors. 
Note that this AC/3 indicator is a promising tool for photometric 
classification of SNe la, provided its scatter remains compara- 
ble to that found for the intermediate redshift sample. 



Sample 


Bands 


Events 


r.m.s 


Average 


nearby 


UBV 


28 


0.122 


0.0008 ± 0.023 


intermediate 




10 


0.033 


0.009 ± 0.010 


high-z 




17 


0.156 


0.039 ± 0.035 



Table 4 Statistics of the 3 samples displayed in Fig.FTTI 



the (r^, i^) pair roughly changes from rest-frame {,B,V) to rest- 
frame (U,B). 

Our cosmological conclusions rely on having a consistent 
distance estimate when using rest- frame BV and UB. This 
property is tested in iGuv et al.l(l2005h . However, it can be tested 
further on the subset of SNLS data having at least three usable 
photometric bands. The test proceeds as follows: 

1 . We fit the three bands at once, and store the stretch and date 
of maximum B light. 



The same exercise can be done without imposing identical 
stretch and date of maximum light on the two fits. Rather than 
testing the light curves model, one then tests for potential biases 
in color estimates (leading to biases in distance estimates). The 
conclusions are the same as with fixed parameters: the sam- 
ples have averages consistent with 0, and the dispersion of the 
central sample increases from 0.033 to 0.036. 
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Fig. 1 1 AU2, difference between rest-frame U peak magnitude 
"predicted" from B and V, and the measured value, as a func- 
tion of redshift. The error bars reflect photometric uncertain- 
ties. The redshift regions have been chosen so that the mea- 
sured bands roughly sample the U BV rest-frame region. The 
differences between average values for the three samples agree 
within statistical uncertainties, indicating that the relation be- 
tween U, B and V brightnesses does not change with redshift. 
Although the nearby and intermediate samples have compara- 
ble photometric resolution, the intermediate sample exhibits a 
far smaller scatter. We attribute this difference to the practical 
difficulties in calibrating U band observations. 



7. Systematic uncertainties 

We present, in this Section, estimates of the systematic un- 
certainties possibly affecting our cosmological parameter mea- 
surements. 

7.1. Photometric calibration and filter band-passes 

We simulated a zero-point shift by varying the magnitudes of 
the light-curve points, one band at a time. Table|5lgives the re- 
sulting shifts in the derived cosmological parameters from the 
calibration eiTors derived in Section RTTI We assume that errors 
in the gMfMiMZM zero-points are independent, and propagate 
these 4 errors quadratically to obtain the total effect on cosmol- 
ogy- 



We rely on the spectrum of one object, Vega (a Lyrae), to 
transform magnitude s into fluxes; th e broadband flux errors for 
Vega are about 1% jHavesI ( 1 19851) and Section B31. To take 
into account the Vega flux and broadband color uncertainties, 
we simulated a flux error linear in wavelength that would offset 
the Vega {B - R) color by 0.01. The impact on D.m is +0.012. 

Uncertainties in the filter bandpasses affect the determina- 
tion of supernovae brightnesses; the first-order effect is from 
errors in the central wavelengths. In the color-color relations 
(Landolt/MegaCam and SDSS/MegaCam - Section l4!2t . we 
were able to detect shifts of 10 A (corresponding roughly to 
a change of 0.01 in the color term). The effect of this shift is 
in fact very small: only the filter has a sizable impact of 
±0.007 on D.M. 

7.2. Light-curve fitting, (U-B) color and k-corrections 

If the light-curve model fails to properly describe the true light- 
curve shape, the result would be a bias in the light-curve pa- 
rameters, and possibly in the cosmological parameters if the 
bias depends on redshift. We have already discussed two possi- 
ble causes of such a bias: the influence of the first measure- 
ment date (Section 15. 2> . and the choice of rest-frame bands 
used to measure brightness and color (Section f6.3> . Both have 
very small effects. However, given only 10 intermediate red- 
shift SNLS events, each with an uncertainty of 0.033, the pre- 
cision with which we can define the average {U - B) color at 
given (B - V) is limited to about 0.01 mag by our sample size. 

Uncertainties in the k-corrections (due to SNe la spectral 
variability at fixed color) contribute directly to the observed 
scatter. The redshift range of the intermediate redshift sam- 
ple of Section l631 corresponds to a rest-frame wavelength span 
of about 400 A, in a region where SNe la spectra are highly 
structured. Since we observe compatible intrinsic dispersions 
for nearby and SNLS events (indeed, slightly lower for SNLS), 
we find no evidence that k-correction uncertainties add signifi- 
cantly to the intrinsic dispersion. 

Nevertheless, since the measured scatter of the intermediate 
redshift sample appears surprisingly small and, since the sam- 
ple size is small, we used a more conservative value of 0.02 
for the light-curve model error, to account for both the errors in 
the colors and from k-corrections. A shift of the {/-band light- 
curve model of 0.02 mag results in a change in Qm of 0.018. 
This is to be added to the statistical uncertainty. 



Band 


zero-point shift 


6Q.U (flat) 




5w (fixedClivi) 


8m 


0.01 


0.000 


-0.02 


0.00 


I'm 


0.01 


0.009 


0.03 


0.02 


'm 


0.01 


-0.014 


0.17 


-0.04 


Zm 


0.03 


0.018 


-0.48 


-0.03 


sum 




0.024 


0.51 


0.05 



Table 5 Influence of a photometric calibration error on the cos- 
mological parameters. 



7.3. U -band variability and evolution of SNe la 

Concerns have been expressed regarding the use of rest-frame 
t/- band fluxes to measu re luminosity distances (e.g. ljhall2002l 
and lNugent et aljr2002l) . motivated by the apparent large vari- 
ability of the t/-band luminosity of SNe la. Such variability 
seems also to be present at intermediate redshifts although 
there seems t o be little ob vious evolution to z = 0.5 of the over- 
afl UV SED dEUis et all in prep.). Note that lGuv et all (l2005l) 
have succeeded in constructing a distance estimator using U 
and B-band data which shows a dispersion of only 0.16 mag 
around the Hubble fine, comparable to that found for distances 
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derived using B- and V-band data. Note also that the quan- 
tity AU3 appears to be independent of redshift, implying that 
if the average luminosity of SNe la evolves with redshift, this 
evolution must pr eserve the UBV rest-frame color relations. 
iLentz et al] ( l2000l) predict a strong dependence of the UV flux 
from the progenitor metallicity (at fixed B - V color), which 
should have been visible if metallicity evolution were indeed 
present. 

7.4. Malmquist bias 

The Malmquist bias may affect the cosmological conclusions 
by altering the average brightness of measured SNe in a red- 
shift dependent way. The mechanism is however not exactly 
straightforward since the reconstructed distance depends on 
stretch and color, and not only on the brightness. We have 
conducted simulations, both of nearby SN searches and of the 
SNLS survey, to investigate the effects on the derivation of cos- 
mological parameters. 

We simulated light-curves of nearby SNe la (0.02 < z < 
0.1) with random explosion date, stretch and color, using the 
observed brighter-slower and brighter-bluer correlations. We 
then simulated a brightness cut at a fixed date. Although the 
number of "detected" events and their average redshift strongly 
depends on the brightness cut, the average distance bias of the 
survivors is found to change by less than 10%, when varying 
both the value and the sharpness of the brightness cut. The bias 
is also essentially independent of the discovery phase, although 
the peak brightness is not. We find a distance modulus bias of 
0.027 (similar in B, V and R), sensitive at the 10% level to the 
unknown details of nearby searches. Note that the redshift de- 
pendence of the distance bias of the nearby sample has no im- 
pact on the cosmological measurements: only the average bias 
matters. 

The crude simulation we conducted applies only to flux 
limited searches, which applies to about half of the sample. We 
compute an average bias value for our nearby sample as the 
simulation result (0.027 mag) times the fraction of events to 
which it applies. Assuming that both factors suffer from an un- 
certainty of 50%, we find an average nearby sample bias value 
of 0.017 ± 0.012 mag. A global increase of all nearby distances 
by 0.017(+0.012) mag increases Qm (flat universe) by 0.019 
(+0.013). 

For the distant SNLS sample, which is flux limited, we sim- 
ulated supernovae at a rate per co-moving volume independent 
of redshift, accounted for the brighter- slower and brighter-bluer 
correlations, and adjusted the position and smoothness of the 
limiting magnitude cut in order to reproduce the redshift and 
peak magnitude distributions. In contrast with nearby SN sim- 
ulations, here we have many observed distributions for a sin- 
gle search, and the key parameters that enter the simulation 
are highly constrained. The best match to SNLS data is shown 
in Figure and Figure shows the expected biases as a 
function of redshift in the shape and color parameters, and for 
our distance estimator The distance modulus bias is about 0.02 
mag at z = 0.8, increasing to 0.05 at z = 1. Correcting for the 
computed bias decreases Q.m (flat Universe) by 0.02. We as- 
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Fig. 12 Distributions of redshifts, peak !m magnitudes (AB), 
stretch factors and colors of SNLS supernovae (black dots) to- 
gether with the distributions obtained with simulated SNe (red 
histograms). 



sumed that the uncertainty in this bias correction is 50% of its 
value. 

To summarize: we find that the differential bias between 
nearby and distant samples almost exactly cancels, and esti- 
mate an overall uncertainty of 0.016 in Qm (flat Universe). 
Since applying the Malmquist bias corrections changes the cos- 
mological results by less than 0. 1 cr, the corrections have not 
been applied. However, in the future, when the SNLS sample 
size increases, modeling and applying the Malmquist bias cor- 
rection will assume a greater importance. The same applies to 
the nearby sample, where having a more controlled and ho- 
mogene ous sample, discovered by a single search (e.g. SN 
Factory, lAlderinget al. 2002) will be essential to reduce the 
associated systematic uncertainty. 

7.5. SNe Ib/c interlopers 

All supernovae used here were spectroscopically identified as 
SN la, but we ha ve labeled the least s ecure identifications as SN 
la* (Section lOFHowell et al.l2"ool . These 15 events are prob- 
able SN la but for this class a small amount of contamination 
by SNe lb or SNe Ic (SNe Ib/c) is possible. We have checked 
that cosmological fits done with or without these events lead to 
the same cosmological conclusions (Section l54l . 

We also looked at estimating the SN Ib/c contamination in 
our sample. SNe Ib/c have an intrinsic luminosity distribution 
which is wider than SNe la (cf. dispersion 0.45 mag fo r SNe la, 
vs. 1 .2 mag for SNe Ib/c; iHomeier 2005: Richardson et alJ 
l2002h . After correcting for the SNe la brighter-slower and 
brighter-bluer correlations, a conservative estimate is that the 
SNe Ib/c scatter around the SNe la Hubble line with a disper- 
sion 3 to 4 times larger than for SNe la. The first clue of SN Ib/c 
contamination would be the presence of objects with large 
residuals around the Hubble line; these contaminants should 
on average be fainter than SNe la at the same redshift. We have 
rejected two objects from the Hubble diagram (Section l5.4> . 
Even if we consider both of these events to be SN Ib/c events, 
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and assume that the dispersion of the SN la distribution about 
the Hubble line is 4 times smaller than for SNe Ib/c, we expect 
on average only 0.5 Ib/c interloper within the fitted sample. 

For these reasons, we estimate the potential bias arising 
from the presence of non la events in our sample to be neg- 
ligible. 

7.6. Gravitational lensing and grey dust 

Gravitational lensing by mass inhomogeneities may affect the 
apparent brightness of our supernovae. With respect to a uni- 
form matter density, most of the events experience a tiny de- 
amplification, and a small minority are amplified (see e.g. 
IhoIz & Wald 1998). 

Whereas the average flux is conserved in the case of weak 
lensing, part of the SN light is lost when strong lensing pro- 
duces multiple images among which some escape detection. 
Multiple images of distant radio sourc es have been syste mat- 
ically searched by the CLASS survey jMvers et al.ll2003l) and 
have proved to be rare: the occuiTence of multiple images sep- 
arated by more than 0.3" and with flux ratio below 10:1 was 
found to be of 1 out of 690 with 1 .44 secondary image on aver- 
age, with inefficiencies due to t he separation and flux ratio cuts 
of 13% and 37% respectively jBrowne et al]l2003h . Multiple 
images with a smaller separation are not resolved in the SNLS, 
and their time delay is much smaller than the typical duration 
of a SN light curve^ so that no flux is lost for such events. 
Hence CLASS results provide us with an upper limit for the 
number of (resolved) strong lensing cases in the SNLS super- 
nova sample, giv en the fact that CLASS sources are globally 
more distant (see lChael2003l) . Assuming (pessimistically) that 
for each strongly lensed SN, we see only one image, the flux 
bias is smaller than 0.3% at z = 1. 

Gravitational lensing also broadens asymmetri- 
cally t he brightness distribu tion of SNe at large red- 
shifts jBer^strom et alJ l2000l) . As a consequence, a cos- 
mologi cal fit using SNe magnitudes (instead of fluxes) is 
biased. Ih oIz & Linder (2004) fo und a dispersion of 0.088 x z 
(note fliat lBergstrom et alJIioOOl find a value of -0.04 at z = 1 
for smooth halo profiles in flat ACDM), which translates into a 
bias of the average magnitude of ~ 0.004 x z- The broadening 
of the brightness distribution also affects the cosmological 
parameters uncertainties. In the cosmological fit, we have 
derived a constant "intrinsic" dispersion which includes the 
average dispersion due to lensing. Neglecting its redshift 
dependence has no significant impact on the accuracy of the 
errors derived for the cosmological parameters. 

In summary, the total effect of lensing on cosmological pa- 
rameters is very small. We find that Qm for a flat ACDM cos- 
mology and the equation of state for a flat universe with BAO 
constraints, are shifted by at most -0.005 and -0.01 respec- 
tively. We therefore did not apply any correction to our results. 

The possibility that SNe la could be dimmed by intergalac- 
tic grey dust (i.e. with weak extinction variation over the optical 

^ Delays are of order of a day for a source at z = 1 and 
a point-like l ens at z = 0.5 for a typical angular separation of 
0.2" tiBergstrom et allZOOOh 



wavelengths) has been suggested bv lAguirrd Jl999albl) as an 
astrophysical alternative to the dark energy hypothesis. Some 
simple dust scenarios without a cosmological constant could be 
excluded bv lRiess et alJ ( l2004l) usin g SNe la data. Studying the 
colors of a large sample of quasars. lOstman & Mortse Ul ll2005h 
were able to set limits on the light absorption length as a func- 
tion of Ry, but these limits can only be translated into an up- 
per bound of supernovae dimming. Conservatively assum ing 
Rv — 12, using the SNOC program llGoobar et al. ] l2002l) . we 
computed an upper limit in the dimming of supernovae which 
translates into a shift of -0.025 in Qm for a ACDM cosmology, 
and a shift of -0.048 in w for a flat cosmology with constant 
equation of state when combined with SDSS BAO results. Note 
that these are upper limits and that a scenario without any in- 
tergalactic dust cannot be excluded. We therefore did not apply 
any correction to our results. 

8. Summary and perspectives 

Table|6lsummarizes the uncertainties affecting our cosmologi- 
cal parameter measurements. The table includes the impact of 
uncertainties in several parameter directions: the Qm direction 
for a flat (Qm, ^^a) (i e- w - -I) cosmology, the Qtot direction 
for a general (Qm, Qa) cosmology, and the w direction at fixed 
Qm for a (Qm, w) cosmology. We also report here the observed 
shifts when the BAO prior is applied to a flat (Qm, w) cosmol- 
ogy- 



Source 




cr(f2,o,) 


cr(vy) 


O-(fiM) 


cr(w) 




(flat) 






(with BAO) 


Zero-points 


0.024 


0.51 


0.05 


0.004 


0.040 


Vega spectrum 


0.012 


0.02 


0.03 


0.003 


0.024 


Filter bandpasses 


0.007 


0.01 


0.02 


0.002 


0.013 


Malmquist bias 


0.016 


0.22 


0.03 


0.004 


0.025 


Sum (sys) 


0.032 


0.55 


0.07 


0.007 


0.054 


Meas. errors 


0.037 


0.52 


0.09 


0.020 


0.087 


U-B color(stat) 


0.020 


0.10 


0.05 


0.003 


0.021 


Sum (stat) 


0.042 


0.53 


0.10 


0.021 


0.090 



Table 6 Summary of uncertainties in the derived cosmological 
parameters. The dominant systematic uncertainty arises from 
the photometric calibration, itself dominated by the im and zm 
band contributions. 



Note that measurement and isolation of systematic errors is a 
major goal of the SNLS. Some of these uncertainties will de- 
crease as more data is acquired and future papers will examine 
a wider range of issues, using our growing dataset. 

Combining Tables|3]and|6j we obtain the following results: 

Qm = 0.263 ± 0.042 (stat) + 0.032 (syst) 
for a flat A cosmology, and 

Qm = 0.271 + 0.021 (stat) + 0.007 (syst) 
w = -1.023 + 0.090 (.?fflf) + 0.054(53/5?) 
w < -0.85 (95% CL) 
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for a flat cosmology with constant equation of state, when 
our constraints are combined with the BAO SDSS results. 
Assuming w > -1 brings our upper limit to -0.83 (at 95% CL). 
Supernovae alone give a marginal constraint: w < -0.5 at 
95% CL. 

These results agree well with previous work s, both from 
SNe la , and also from other sources. For example lSeliak et al] 
( l2005h finds very similar results combining CMB, LSS and 
Lyo- constraints. The dominant systematic errors arise from the 
nearby sample and from the photometric calibration of the zm 
band; both will be improved in the future. The multi-band light- 
curves allow us to study color relations as a function of redshift; 
these data are expected to be sensitive indicators of evolution. 
We observed a surprisingly narrow correlation between (U-B) 
and (B-V) (using the Af/j indicator), indicating that the disper- 
sion in U -band properties is well correlated with measurements 
in redder bands. 

From the first year of SNLS data, we placed 71 distant 
events on the Hubble Diagram, with 10 more from the same 
period to be added later (Our full first year statistics would 
have been around 100 SNe la with spectroscopic confirmation 
had we not lost Feb 2004 to an instrument failure.) Our 
time sampling, filter coverage, and image quality have now 
significantly improved since early 2004, and we now regularly 
acquire 2-3 times as much data m zm- A precise photometric 
calibration is essential, and we are now working with the 
CFHTLS community in refining the photometric calibration 
of the MegaCam instrument. We have embarked on the 
process of calibrating tertiary standards in our fields, from 
Sloan secondary and primary standards. This will allow us to 
cross-check the Vega/Landolt zero-points, and more accurately 
calibrate zm band observations. 

After only two years of operation, the SNLS has already 
demonstrated its advantages over all previous ground-based 
supernova surveys. The "rolling search" technique is robust 
to weather and instrument-related problems, and the technical 
characteristics of the survey are now well understood. The av- 
erage rate of spectroscopically-confirmed SNe la is currently 
about 10 per lunation and continues to increase. Up until July 
2005, the SNLS sample includes more than 200 spectroscopi- 
cally identified SNe la, most with excellent photometric tem- 
poral and filter coverage. An extrapolation of the current rate 
to the end of the survey indicates that we should reach our goal 
of building a Hubble diagram with about 700 spectroscopically 
identified well-measured SN la events. The SNLS already has 
the largest-ever sample of high-z SNe discovered by a single 
telescope, and will eventually produce a homogeneous, high- 
quaUty sample that is an order of magnitude larger still. 

High statistical accuracy benefits the control of systematics. 
With our unmatched SN statistics, by year 5 we will be able to 
populate each ~ 0. 1 redshift bin with ~ 100 SNe la, thus filling 
the brightness, decline-rate, and color 3-dimensional parameter 
space. This will enable us to detect possible drifts in "SNe la 
demographics", and control Malmquist bias. Moreover, the 
rolUng-search observing mode produces many events at low 
to intermediate redshift with superb photometric accuracy, be- 
cause integration times are tailored for the faintest objects. 



These relatively bright events permit demanding internal con- 
sistency tests, and may lead to improvements in distance esti- 
mation. 
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Fig. 13 Stretch, color and Hubble diagram residuals as a func- 
tion of redshift for SNLS supernovae (gray dots). The black 
points correspond to average values in redshift bins. The red 
soUd (dashed) Unes represent the average (one standard devi- 
ation) values obtained with SNe simulations as described in 
Section lT4l At large redshifts, since only bright SNe are iden- 
tified, the average stretch factor is larger and the average color 
bluer. The average distance modulus is less affected by the se- 
lection (see text for details). 
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Table 7. Transients from the first year sample identified as SNIa or SNIa* 



Name 


RA(2000) 


Dec(2000) 


Obs. date" 


Tel./Inst.'' 


Spectral Type"^ 




z source 


SNLS- 


-03Dlau 


02:24:10.392 


-04:02:14.93 


2907 


Keck/LRIS 


SNIa 


0.504 


gal 


SNLS- 


-03Dlaw 


02:24:14.786 


-04:31:01.61 


2907 


Keck/LRIS 


SNIa 


0.582 


gal 


SNLS- 


-03Dlax 


02:24:23.338 


-04:43:14.28 


2913 


Gem/GMOS 


SNIa 


0.496 


gal 


SNLS- 


-03Dlbf 


02:24:02.375 


-04:55:57.27 


2909 


VLT/FORSl 


SNIa* 


0.703 


gal 


SNLS- 


-03Dlbk 


02:26:27.412 


-04:32:11.95 


2912 


Gem/GMOS 


SNIa 


0.865 


gal 


SNLS- 


-03Dlbp 


02:26:37.714 


-04:50:19.55 


2910 


VLT/FORSl 


SNIa* 


0.346 


gal 


SNLS- 


-03Dlcm 


02:24:55.294 


-04:23:03.61 


2940 


Gem/GMOS 


SNIa 


0.87 


SN 


SNLS- 


-03D1CO 


02:26:16.252 


-04:56:05.65 


2966 


Keck/LRIS 


SNIa 


0.679 


gal 


SNLS- 


-03Dldj 


02:26:19.087 


-04:07:08.89 


2964 


Keck/LRIS 


SNIa 


0.39 


SN 


SNLS- 


-03Dldt 


02:26:31.200 


-04:03:08.51 


2974 


VLT/FORSl 


SNIa 


0.612 


gal 


SNLS- 


-03Dlew 


02:24:14.079 


-04:39:56.93 


2995 


Gem/GMOS 


SNIa 


0.868 


gal 


SNLS- 


-03Dlfb 


02:27:12.875 


-04:07:16.44 


2641 


VLT/FORSl 


SNIa 


0.498 


gal 


SNLS- 


-03Dlfc 


02:25:43.625 


-04:08:38.93 


2641 


VLT/FORSl 


SNIa 


0.331 


gal 


SNLS- 


-03Dlfl 


02:25:58.329 


-04:07:44.17 


2641 


VLT/FORSl 


SNIa 


0.688 


gal 


SNLS- 


-03Dlfq 


02:26:55.683 


-04:18:08.10 


2998 


Gem/GMOS 


SNIa 


0.80 


SN 


SNLS- 


-03Dlgt 


02:24:56.027 


-04:07:37.11 


2641 


VLT/FORSl 


SNIa 


0.55 


SN 


SNLS- 


-03D3af 


14:21:14.955 


+52:32:15.68 


2737 


Keck/LRIS 


SNIa 


0.532 


gal 


SNLS- 


-03D3aw 


14:20:53.534 


+52:36:21.04 


2767 


Keck/LRIS 


SNIa 


0.449 


gal 


SNLS- 


-03D3ay 


14:17:58.448 


+52:28:57.63 


2766 


Keck/LRIS 


SNIa 


0.371 


gal 


SNLS- 


-03D3ba 


14:16:33.465 


+52:20:32.02 


2766 


Keck/LRIS 


SNIa 


0.291 


gal 


SNLS- 


-03D3bb 


14:16:18.920 


+52:14:53.66 


2766 


Keck/LRIS 


SNIa 


0.244 


gal 


SNLS- 


-03D3bh 


14:21:35.894 


+52:31:37.86 


2766 


Keck/LRIS 


SNIa 


0.249 


gal 


SNLS- 


-03D3bl 


14:19:55.844 


+53:05:50.91 


2792 


Keck/LRIS 


SNIa 


0.355 


gal 


SNLS- 


-03D3CC 


14:19:45.192 


+52:32:25.76 


2793 


Keck/LRIS 


SNIa 


0.463 


gal 


SNLS- 


-03D3cd 


14:18:39.963 


+52:36:44.22 


2792 


Keck/LRIS 


SNIa 


0.461 


gal 


SNLS- 


-03D4ag 


22:14:45.806 


-17:44:22.95 


2824 


Keck/LRIS 


SNIa 


0.285 


gal 


SNLS- 


-03D4at 


22:14:24.023 


-17:46:36.03 


2826 


VLT/FORSl 


SNIa 


0.633 


gal 


SNLS- 


-03D4au 


22:16:09.917 


-18:04:39.19 


2826 


VLT/FORSl 


SNIa* 


0.468 


gal 


SNLS- 


-03D4bc 


22:15:28.143 


-17:49:48.66 


2826 


VLT/FORSl 


SNIa 


0.572 


gal 


SNLS- 


-03D4cj 


22:16:06.658 


-17:42:16.83 


2879 


Keck/LRIS 


SNIa 


0.27 


SN 


SNLS- 


-03D4cn 


22:16:34.600 


-17:16:13.55 


2879 


Gem/GMOS 


SNIa 


0.818 


gal 


SNLS- 


-03D4CX 


22:14:33.754 


-17:35:15.35 


2885 


VLT/FORSl 


SNIa 


0.95 


SN 


SNLS- 


-03D4cy 


22:13:40.441 


-17:40:54.12 


2909 


Gem/GMOS 


SNIa 


0.927 


gal 


SNLS- 


-03D4CZ 


22:16:41.845 


-17:55:34.40 


2910 


Gem/GMOS 


SNIa 


0.695 


gal 


SNLS- 


-03D4dh 


22:17:31.040 


-17:37:46.98 


2906 


Keck/LRIS 


SNIa 


0.627 


gal 


SNLS- 


-03D4di 


22:14:10.249 


-17:30:24.18 


2885 


VLT/FORSl 


SNIa 


0.905 


gal 


SNLS- 


-03D4dy 


22:14:50.513 


-17:57:23.24 


2912 


VLT/FORSl 


SNIa 


0.60 


SN 


SNLS- 


-03D4fd 


22:16:14.462 


-17:23:44.33 


2937 


Gem/GMOS 


SNIa 


0.791 


gal 


SNLS- 


-03D4gf 


22:14:22.907 


-17:44:02.49 


2641 


VLT/FORSl 


SNIa* 


0.58 


SN 


SNLS- 


-03D4gg 


22:16:40.185 


-18:09:51.82 


2641 


VLT/FORSl 


SNIa 


0.592 


gal 


SNLS- 


-03D4gl 


22:14:44.183 


-17:31:44.36 


2966 


Keck/LRIS 


SNIa 


0.571 


gal 


SNLS- 


-04Dlag 


02:24:41.125 


-04:17:19.66 


2641 


VLT/FORSl 


SNIa 


0.557 


gal 


SNLS- 


-04Dlaj 


02:25:53.982 


-04:59:40.50 


2641 


VLT/FORSl 


SNIa* 


0.72 


SN 


SNLS- 


-04Dlak 


02:27:33.399 


-04:19:38.73 


2641 


VLT/FORSl 


SNIa* 


0.526 


gal 


SNLS- 


-04D2ac 


10:00:18.923 


+02:41:21.63 


2641 


VLT/FORSl 


SNIa* 


0.348 


gal 


SNLS- 


-04D2ae 


10:01:52.361 


+02:13:21.27 


3026 


Gem/GMOS 


SNIa 


0.843 


gal 


SNLS- 


-04D2al 


10:01:52.482 


+02:09:51.25 


2641 


VLT/FORSl 


SNIa 


0.84 


SN 


SNLS- 


-04D2an 


10:00:52.332 


+02:02:28.73 


2641 


VLT/FORSl 


SNIa 


0.62 


SN 


SNLS- 


-04D2bt 


09:59:32.725 


+02:14:53.07 


3085 


VLT/FORSl 


SNIa 


0.220 


gal 
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Name 


RA(2000) 


Dec(2000) 


Obs. date^ 


Tel./Inst.'' 


Spectral Type'^ 




z source 


SNLS- 


-04D2ca 


10:01:20.514 


+02:20:21.76 


2641 


VLT/FORSl 


SNIa 


0.83 


SN 


SNLS- 


-04D2cf 


10:01:56.110 


+01:52:46.40 


3086 


VLT/FORSl 


SNIa 


0.369 


gal 


SNLS- 


-04D2CW 


10:01:22.787 


+02:11:55.31 


3085 


VLT/FORSl 


SNIa 


0.568 


gal 


SNLS- 


-04D2fp 


09:59:28.162 


+02:19:15.58 


2641 


VLT/FORSl 


SNIa 


0.415 


gal 


SNLS- 


-04D2fs 


10:00:22.110 


+01:45:55.70 


2641 


VLT/FORSl 


SNIa 


0.357 


gal 


SNLS- 


-04D2gb 


10:02:22.676 


+01:53:39.34 


3117 


Keck/LRIS 


SNIa 


0.43 


SN 


SNLS- 


-04D2gc 


10:01:39.281 


+01:52:59.36 


3118 


Keck/LRIS 


SNIa 


0.521 


gal 


SNLS- 


-04D2gp 


09:59:20.400 


+02:30:31.88 


3116 


VLT/FORSl 


SNIa 


0.71 


SN 


SNLS- 


-04D2iu 


10:01:13.221 


+02:24:53.91 


3139 


VLT/FORSl 


SNIa* 


0.69 


SN 


SNLS- 


-04D2ja 


09:58:48.519 


+01:46:18.64 


3139 


VLT/FORSl 


SNIa* 


0.74 


SN 


SNLS- 


-04D3bf 


14:17:45.096 


+52:28:04.31 


3054 


Gem/GMOS 


SNIa 


0.156 


gal 


SNLS- 


-04D3CO 


14:17:50.030 


+52:57:48.95 


3117 


Keck/LRIS 


SNIa 


0.62 


SN 


SNLS- 


-04D3cp 


14:20:23.954 


+52:49:15.45 


3119 


Keck/LRIS 


SNIa 


0.83 


SN 


SNLS- 


-04D3cy 


14:18:12.452 


+52:39:30.40 


3115 


Keck/DMOS 


SNIa 


0.643 


gal 


SNLS- 


-04D3dd 


14:17:48.411 


+52:28:14.57 


3122 


Gem/GMOS 


SNIa 


1.01 


SN 


SNLS- 


-04D3df 


14:18:10.042 


+52:16:39.85 


3117 


Keck/LRIS 


SNIa 


0.47 


SN 


SNLS- 


-04D3do 


14:17:46.113 


+52:16:03.36 


3117 


Keck/LRIS 


SNIa 


0.61 


SN 


SNLS- 


-04D3ez 


14:19:07.894 


+53:04:19.17 


3118 


Keck/LRIS 


SNIa 


0.263 


gal 


SNLS- 


-04D3fk 


14:18:26.198 


+52:31:42.74 


3118 


Keck/LRIS 


SNIa 


0.358 


gal 


SNLS- 


-04D3fq 


14:16:57.902 


+52:22:46.46 


3123 


Gem/GMOS 


SNIa 


0.73 


SN 


SNLS- 


-04D3gt 


14:22:32.611 


+52:38:49.30 


3149 


Keck/LRIS 


SNIa* 


0.451 


gal 


SNLS- 


-04D3gx 


14:20:13.666 


+52:16:58.33 


3147 


Gem/GMOS 


SNIa* 


0.91 


SN 


SNLS- 


-04D3hn 


14:22:06.908 


+52:13:43.00 


3148 


Gem/GMOS 


SNIa 


0.552 


gal 


SNLS- 


-04D3is 


14:16:51.968 


+52:48:45.70 


3149 


Keck/LRIS 


SNIa* 


0.71 


SN 


SNLS- 


-04D3ki 


14:19:34.598 


+52:17:32.61 


3149 


Keck/LRIS 


SNIa* 


0.930 


gal 


SNLS- 


-04D3ki- 


14:16:35.943 


+52:28:44.02 


3173 


Gem/GMOS 


SNIa 


0.337 


gal 


SNLS- 


-04D3ks 


14:22:33.479 


+52:11:07.44 


3149 


Keck/LRIS 


SNIa* 


0.752 


gal 


SNLS- 


-04D31p 


14:19:50.911 


+52:30:11.88 


3153 


Gem/GMOS 


SNIa* 


0.983 


gal 


SNLS- 


-04D31U 


14:21:08.009 


+52:58:29.74 


3180 


Gem/GMOS 


SNIa 


0.822 


gal 


SNLS- 


-04D3mk 


14:19:25.768 


+53:09:49.48 


3176 


Gem/GMOS 


SNIa 


0.813 


gal 


SNLS- 


-04D3ml 


14:16:39.095 


+53:05:35.89 


3177 


Gem/GMOS 


SNIa 


0.95 


SN 


SNLS- 


-04D3nc 


14:16:18.224 


+52:16:26.09 


3200 


Gem/GMOS 


SNIa* 


0.817 


gal 


SNLS- 


-04D3nh 


14:22:26.729 


+52:20:00.92 


3180 


Gem/GMOS 


SNIa 


0.340 


gal 


SNLS- 


-04D3nq 


14:20:19.193 


+53:09:15.90 


3201 


Gem/GMOS 


SNIa 


0.22 


SN 


SNLS- 


-U4U3nr 


14:22:38.526 


+52:38:55.89 


3202 


Gem/GMOS 


SNIa* 


0.96 


SN 


SNLS- 


-04D3ny 


14:18:56.332 


+52:11:15.06 


3197 


Gem/GMOS 


SNIa 


0.81 


SN 


SNLS- 


-04D3oe 


14:19:39.381 


+52:33:14.21 


3198 


Gem/GMOS 


SNIa 


0.756 


gal 


SNLS- 


-04D4an 


22:15:57.119 


-17:41:43.93 


3200 


VLT/FORSl 


SNIa 


0.613 


gal 


SNLS- 


-04D4bk 


22:15:07.681 


-18:03:36.79 


3200 


VLT/FORSl 


SNIa* 


0.84 


SN 


SNLS- 


-04D4bq 


22:14:49.391 


-17:49:39.37 


3203 


VLT/FORSl 


SNIa 


0.55 


SN 


SNLS- 


-04D4dm 


22:15:25.470 


-17:14:42.71 


3206 


Gem/GMOS 


SNIa 


0.811 


gal 


SNLS- 


-04D4dw 


22:16:44.667 


-17:50:02.38 


3206 


VLT/FORSl 


SNIa 


0.96 


SN 



"Date of spectrocopic observations (ID 2450000+). 

''Telescope and instrument with which the spectrum was acquired. 

■^See Sect. l2.2l for detinitions. 

■'SN spectrum (SN) or host galaxy spectrum (gal). Sz ~ 0.01 when from SN spectrum, ~ 0.001 when from host galaxy 
spectrum 
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Table 8. Nearby Type la supernovae 



Name 


z 


Bands 






'b 








s 




c 
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U, 
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u. 
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1 
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1^ 
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n 

u. 
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lo. 


.jyZ. 




1 
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n 

u. 
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10 


7A1 
.Z^l 
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n 

u. 
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yrVyO) 
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1 Q 

ly 
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.Lyy 
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U. 
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-U.U^o ± 
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IS 
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1 J 
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u.uu / 


1 

1, 
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u.uu / 
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: U.UUo 


lA 
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u. 
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1 7 
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U.UIO 


n 

U. 
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U.UIO 
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lA 
JO. 
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u. 
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.Uj / 
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1 7 
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U.Uj J 


n 

U. 


7SA 
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U.UIO 
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-U-Ul^ ± 
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lA 
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.Z /O ± 


n 

u. 


OSQ 

.UJ7 
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U.Ul o 
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U.UiZ 


U. 
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u.uuo 
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U.Ujj ± 
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j4. 
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u. 
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1 o 
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U. 
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n nAi a. 

-U.U^J ± 
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17 
J / . 
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n 

u. 
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.U^l 






U.Uoo 


R\A 
D V 


1 Q 

Ly 


.jyo 




U.U / J 


n 

U. 


A^n 




n 07 Qd 

u.uzy 


017 -1- 
U.UjZ ± 


n ni7d 

U.Uo / 


IS 
JO. 


171 -1- 
. IZl ± 


u. 


OAA 






U.UOj 


R\A 
D V 


lo 


1 77 




U.U'^-l 


1 

1, 


nni 

.UUl 




01 sd 
u.uio 


n niA a. 


n ni Qd 

: U.UI7 


17 

J / . 


SAO -1- 


u. 


OAA 




1 QQ9P 
ly yZr 


u.uzo 


RV 
D V 


10. 


.Uj / 




U.uio 


1 

1 . 


1 1Q 

. 1 J7 




U.UoH- 


n nos -^ 

-U.UUJ ± 


- n n 1 Rd 
. u.uio 


l'^ 

J J. 


SAS -1- 
.JDJ It 


n 

U. 


1 Al 
. IH- 1 


yrlyvi) 




U.UjU 


R\7 

r> V 


1 7 

1 / , 


7QQ 

. / yy 




n ni A 


U. 


Q 1 '^ 

,7 1 J 




001s 

U.Ul 


OQA -1- 
U.U7O ± 


n ni 7 

U.Ul / 


lA 
JO. 


S77 -1- 
.OZ / ± 


n 

u. 


OAO 
.UOU 






n n7 1 

U.U / i 


R\7 
D V 


1 o 


.Oil 






U. 


QSS 
.700 




077d 
U.UZZ 


OA1 -1- 
U.U'+l ± 


n n7Ad 
u.uzo 


17 

J / . 


AOA -1- 
.OU^ ± 


n 

u. 


OAS 
.U^O 




Lyyjrl 


U.UZj 


RV 
D V 


10. 


7'^'^ 
. / jD 




n ni 7 

U.Ul / 


n 

u. 


AQQ 

.077 




017 
U.UIZ 


7S0 -1- 
U.ZjU ± 


n ni s 

U.Ul J 


IS 

J J. 


1 Q7 -1- 
. 1 7Z It 


n 

u. 


0Q7 
.U7Z 


('T4QA AOA"! 
(^rl70,/\U'4J 


1 77 


U.UJ J 


RV 
D V 


1 7 

1 / , 


Al A 
.D It- 




n ni 1 

U.Ul 1 


n 

u. 


.7U1 




010 
U.UIU 


n ni A -^ 

-U-UIH- ± 


■ n ni 1 

. U.Ul 1 


lA 

JO. 


IQA -V- 
. 1 y^ ± 


n 

u. 


0A7 




i yyHNl 




R\A 
D V 


10 


.ZUD 




n c\A 1 

U.U'^-l 


U. 


S^A 




01 Qd 
U.UI7 


OAO -1- 
U.U'+U ± 


n n77d 

U.UZZ 
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JJ. 


77S -1- 
.ZZO ± 


u. 
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RV 
D V 
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07A 
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U.UIO ± 


n ni 7 

U.Ul / 


IJ. 


071 -1- 
.U / 1 It 


n 

u. 


1 AA 




1 QQ^or' 
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u.u^y 


R\7 
D V 


1 7 

1 / , 


07 A 
.UZO 




n nnQ 
u.uu 7 


1 

1 , 


nA7 
.u^z 
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U.Ul J 


01 -1- 
U.UIU ± 


n ni n 

U.UIU 


lA 
JO. 


ISl -1- 
.JoJ ± 


U. 


OS 1 

.UJ 1 


t^K77,/\U^ ) 


1 QQ^Krl 


u.uio 


R\7 
D V 


1 JZ 

I J 


7AA 
.Z^O 




n nnQ 
u.uuv 


U. 


QQ7 

.yyL 
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u.uu 7 


701 -1- 
U.Z7J ± 
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U.UUo 


lA 
J4. 


OSl -1- 
.UoJ ± 


n 

u. 


1 18 

. 1 JO 


tK77,/\U^^ 


1 7 7Dau 
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RV 
D V 


I Q 

ly. 
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n 

u. 
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-U.U / ± 


■ n ni sd 
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u. 
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10 
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n ni n 
u.uiu 


U. 
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LyyODO 
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u.uuo 
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n no7 
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lA 
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u. 
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LyyVDv 


n ni 7 

U.Ul / 
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I J. 


. JoU 
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U.Ul 7 


n 

u. 
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77'^ -1- 
U.ZZJ It 


n noQd 

U.UU7 


Id 


1 1 Q -1- 
. J 1 7 it 


n 

u. 


1 11 

. 1 J J 






U.UJU 
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r> V 


10 


A'?A 
.OjO 




n n7Qd 

U.UZ7 


1 
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nA-^ 
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IS 
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S77 -1- 
.OZZ ± 


n 

u. 


7 1 

.ZIU 




Lyy 1 Ug 


U.UjU 


T TR\/ 
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10 
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U. 
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U.UZ4 
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U.UUj ± 
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U.UIU 


IS 

J J. 


QQA -1- 
.77^ ± 
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u. 
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.UoU 
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I QQ7V 
Lyy 1 1 


n ni 7 
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TTRV 
UD V 


1 C 

1 J. 


.ZoH- 




U.uzu 


n 

u. 
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.7 10 
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U.UZH- 


008 -1- 

u.uuo ± 
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IJ. 
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u. 
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.868 


± 


0.024^^ 
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.113 + 


0, 
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0.023 


UBV 


16 


.544 


± 


0.007 





.732 


± 


0.006 


0.190 ± 


0.006 


35. 


.146 ± 


0, 
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(J02) 


2000dk 


0.016 


UBV 


15 
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± 
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± 
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.129 ± 


0, 
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15 
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± 
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± 
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0, 
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BV 
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1, 
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0. 
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15 


.271 


± 
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0, 
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± 
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.118 ± 


0, 
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UBV 
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1, 


.004 
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0.007 


34. 
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0, 
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''CMB-centric redshift. 

''Computed with Hq = 70 km s^'Mpc"' . Uncertainty only accounts for photometric uncertainties. 

'^Photo metry Referen ces : H96: iHamuv et all <199d). R99: iRiess et all <1999h . KOI : iKrisciunas et all j200lh . 
J02: |jh3 j2002h A04: lAltavilla et all <2004ir K04a: IKrisciunas et all j2004 j) . K04b: IKrisciunas et^llOOm . 
S02 lStrolgeretalH2002h . 

''First photometric measurement after S-band maximum, see discussion in section lsTI 
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Table 9. SNLS Type la supernovae 



Name z " Bands 



SNLS 


-03Dlau 


SNLS 


-03Dlaw 


SNLS 


-03Dlax 


SNLS 


-03Dlbp 


SNLS 


-03Dlcm 


SNLS 


-03D1CO 


SNLS 


-03Dlew 


SNLS 


-03Dlfc 


SNLS 


-03Dlfl 


SNLS 


-03Dlfq 


SNLS 


-03Dlgt 


SNLS 


-03D3af 


SNLS 


-03D3aw 


SNLS 


-03D3ay 


SNLS 


-03D3ba 


SNLS 


-03D3bh 


SNLS 


-03D3CC 


SNLS 


-03D3cd 


SNLS 


-03D4ag 


SNLS 


-03D4at 


SNLS 


-03D4au'' 


SNLS 


-03D4bc'' 


SNLS 


-03D4cn 


SNLS 


-03D4CX 


SNLS 


-03D4cy 


SNLS 


-03D4CZ 


SNLS 


-03D4dh 


SNLS 


-03D4di 


SNLS 


-03D4dy 


SNLS 


-03D4fd 


SNLS 


-03D4gf 


SNLS 


-03D4gg 


SNLS 


-03D4gl 


SNLS 


-04Dlag 


SNLS 


-04Dlaj 


SNLS 


-04Dlak 


SNLS 


-04D2cf 


SNLS 


-04D2fp 


SNLS 


-04D2fs 


SNLS 


-04D2gb 


SNLS 


-04D2gc 


SNLS 


-04D2gp 


SNLS 


-04D2iu 


SNLS 


-04D2ja 


SNLS 


-04D3CO 


SNLS 


-04D3cp 


SNLS 


-04D3cy 


SNLS 


-04D3dd 


SNLS 


-04D3df 



0.504 


riz 


0.582 


riz 


0.496 


riz 


0.346 


riz 


0.870 




0.679 




0.868 




0.331 




0.688 




0.800 




0.548 




0.532 




0.449 




0.371 




0.291 




0.249 




0.463 




0.461 




0.285 




0.633 




0.468 


pH 7 


0.572 


pH 7 


0.818 


pH 7 


0.949 


pri7 


0.927 


pri7 


0.695 


pn 7 


0.627 


pn 7 


0.905 


pri7 


0.604 


pri7 


0.791 


pri7 


0.581 


pn7 


0.592 


pri7 


0.571 


pH 7 


0.557 


pri7 


0.721 


pri 7 


0.526 


pri7 


0.369 


pri7 


0.415 


pH 7 


0.357 


pn 7 




griz 


0.521 


griz 


0.707 


griz 


0.691 


griz 


0.741 


griz 


0.620 


griz 


0.830 


griz 


0.643 


griz 


1.010 


griz 


0.470 


griz 



22 


978 


+ 





010 


23 


599 


+ 





020 


22 


957 


+ 





oil 


22 


465 


± 





014 


24 


469 


± 





066 


24 


094 


± 





033 


24 


359 


± 





078 


21 


800 


+ 





005 


23 


629 


+ 





015 


24 


519 


± 





030 


24 


119 


+ 





048 


23 


470 


+ 





027 


22 


552 


+ 





016 


22 


201 


+ 





016 


22 


049 


+ 





034 


21 


132 


+ 





018 


22 


558 


+ 





111 


22 


562 


± 





017 


21 


237 


+ 





005 


23 


746 


+ 





020 


23 


856 


+ 





020 


24 


596 


± 





061 


24 


652 


± 





051 


24 


504 


+ 





083 


24 


718 


+ 





109 


24 


019 


+ 





036 


23 


389 


+ 





Oil 


24 


288 


+ 





068 


23 


313 


+ 





010 


24 


212 


+ 





025 


23 


351 


+ 





013 


23 


403 


+ 





024 


23 


269 


+ 





026 


23 


003 


+ 





Oil 


23 


901 


+ 





030 


23 


631 


+ 





028 


22 


340 


+ 





007 


22 


528 


+ 





010 


22 


422 


+ 





008 


22 


796 


± 





018 


23 


321 


+ 





014 


24 


151 


+ 





047 


24 


258 


+ 





048 


24 


098 


+ 





045 


23 


781 


+ 





022 


24 


235 


± 





063 


23 


798 


+ 





021 


25 


120 


+ 





192 


23 


465 


+ 





010 



stretch'' 



1 


124 


+ 





019 


1 


002 


+ 





024 





899 


± 





010 





880 


± 





007 


1 


173 


± 





061 





975 


+ 





032 


1 


028 


+ 





040 





937 


+ 





005 





999 


± 





024 





806 


± 





052 





856 


± 





042 





907 


+ 





023 





955 


± 





013 





968 


± 





010 


1 


036 


+ 





021 





993 


+ 





008 


1 


074 


± 





031 


1 


131 


+ 





034 


1 


059 


+ 





005 





989 


± 





029 


1 


000 


± 





030 





774 


± 





048 





743 


± 





059 





882 


± 





019 


1 


031 


± 





052 





729 


± 





024 


1 


061 


± 





013 


1 


103 


+ 





041 


1 


056 


+ 





001 





919 


+ 





033 


1 


009 


+ 





026 





966 


± 





049 





957 


+ 





033 





944 


+ 





013 


1 


074 


+ 





067 





824 


± 





021 





895 


± 





003 





964 


+ 





010 





942 


± 





009 





777 


± 





013 


1 


065 


± 





024 





801 


± 





002 





800 


± 





035 





945 


± 





036 





895 


± 





017 


1 


110 


+ 





035 





963 


± 





016 


1 


088 


± 





074 





730 


± 





010 



color'' 



0.030 ± 0.018 
0.018 ±0.030 
-0.044 ± 0.021 
0.143 ±0.017 
-0.035 ±0.143 
-0.021 ± 0.047 
-0.102 ±0.169 
0.042 ± 0.004 
-0.070 ± 0.021 
0.027 ± 0.030 
0.244 ± 0.050 
0.029 ± 0.031 
-0.048 ± 0.019 
-0.018 ±0.014 
0.263 ± 0.015 
-0.090 ± 0.013 
-0.070 ± 0.050 
0.025 ±0.011 
-0.061 ± 0.004 
-0.060 ± 0.030 
0.291 ± 0.034 
0.025 ± 0.078 
0.023 ± 0.158 
0.080 ±0.124 
-0.305 ±0.174 
-0.069 ± 0.043 
0.028 ± 0.016 
0.029 ± 0.120 
0.122 ±0.015 
0.028 ± 0.044 
-0.056 ± 0.024 
0.062 ± 0.035 
0.030 ± 0.028 
-0.182 ±0.017 
0.072 ± 0.038 
0.018 ± 0.033 
0.002 ± 0.010 
0.006 ± 0.015 
0.128 ± 0.008 
-0.008 ± 0.025 
0.185 ± 0.022 
-0.052 ± 0.060 
0.074 ± 0.056 
-0.067 ± 0.043 
-0.064 ± 0.030 
-0.448 ±0.180 
0.017 ± 0.029 
-0.071 ±0.205 
0.060 ± 0.017 



42 


429 


+ 





039 


42 


881 


+ 





054 


42 


180 


+ 





038 


41 


367 







021 


44 


095 


+ 





301 


43 


398 


+ 





088 


43 


871 


+ 





344 


40 


946 


+ 





013 


43 


046 


+ 





049 


43 


490 


+ 





090 


42 


825 


+ 





080 


42 


592 


+ 





083 


41 


866 


+ 





044 


41 


488 


+ 





030 


40 


999 


+ 





033 


40 


571 


+ 





020 


42 


089 


+ 





034 


42 


031 


+ 





058 


40 


731 


+ 





015 


43 


133 


+ 





064 


42 


708 


+ 





069 


43 


521 


+ 





135 


43 


532 


+ 





304 


43 


507 


+ 





272 


44 


553 


+ 





380 


43 


023 


+ 





086 


42 


746 


+ 





035 


43 


708 


+ 





258 


42 


515 


+ 





029 


43 


353 


+ 





076 


42 


761 


+ 





047 


42 


562 


+ 





090 


42 


465 


+ 





070 


42 


511 


+ 





029 


43 


209 


+ 





106 


42 


644 


+ 





055 


41 


485 


+ 





016 


41 


772 


+ 





027 


41 


441 


+ 





018 


4.1 






u 


KJJO 


42 


439 


± 





054 


43 


237 


± 





129 


43 


144 


± 





136 


43 


427 


± 





117 


43 


030 


± 





060 


44 


414 


± 





347 


43 


023 


± 





059 


44 


673 


± 





533 


42 


268 


± 





032 



24 
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Table 9 (cont'd) 



Name 




z " 


Bands 




m* 
'"b 








stretch'' 


color'' 










SNLS- 


-04D3do 


0, 


.610 


griz 


23, 


.574 ± 





.014 


0. 


.862 


± 0.013 


-0.079 ± 


: 0.019 


42, 


.19b ± 





.039 


SNLS- 


-04D3ez 





.263 


griz 


21, 


.678 ± 


0, 


.004 


0. 


.895 


± 0.006 


0.091 ± 


0.003 


40, 


.682 ± 





.013 


SNLS- 


-04D3fk 





.358 


griz 


22, 


.532 ± 


0, 


.005 


0. 


.913 


+ 0.005 


0.149 ± 


0.006 


41, 


.474 ± 





.013 


SNLS- 


-04D3fq 





.730 


griz 


24, 


.128 + 





.026 


0. 


.900 


± 0.014 


-0.002 ± 


: 0.037 


43, 


.287 ± 





.075 


SNLS- 


-04D3gt 





.451 


griz 


23, 


.235 ± 





.010 


0. 


.953 


± 0.010 


0.276 ± 


0.016 


42, 


.038 ± 





.030 


SNLS- 


-04D3gx 





.910 


griz 


24, 


.708 ± 





.094 


0. 


.952 


± 0.047 


-0.202 ± 


: 0.163 


44, 


.259 ± 


0, 


.346 


SNLS- 


-04D3hn 





.552 


griz 


23, 


.475 ± 





Oil 


0. 


.898 


±0.011 


0.106 ± 


0.017 


42, 


.461 ± 





.035 


SNLS- 


-04D3is 





.710 


griz 


24, 


.256 ± 





.027 


0. 


.972 


± 0.002 


0.220 ± 


0.038 


43, 


.176 ± 





.077 


SNLS- 


-04D3ki 





.930 


griz 


24, 


.871 ± 





.126 


0. 


.901 


± 0.039 


-0.256 ± 


: 0.194 


44, 


.430 ± 





.430 


SNLS- 


-04D3kr 





.337 


griz 


21, 


.967 ± 





.003 


1. 


.064 


± 0.004 


0.072 ± 


0.003 


41, 


.259 ± 





.010 


SNLS- 


-04D3ks 





.752 


griz 


23, 


.882 ± 





.035 


1. 


.013 


± 0.037 


0.026 ± 


0.043 


43, 


.170 ± 





.090 


SNLS- 


-04D31p 





.983 


griz 


24, 


.925 ± 





.168 


0. 


.831 


± 0.049 


0.022 ± 


0.211 


43, 


.941 ± 





.496 


SNLS- 


-04D31U 





.822 


griz 


24, 


.342 ± 





.040 


0. 


.950 


± 0.028 


0.019 ± 


0.116 


43, 


.544 ± 





.218 


SNLS- 


-04D3ml 





.950 


griz 


24, 


.552 ± 


0, 


.082 


1, 


.182 


±0.015 


0.117 ± 


0.122 


43, 


.954 ± 





.268 


SNLS- 


-04D3nc 





.817 


griz 


24, 


.271 ± 





.048 


1. 


111 


± 0.064 


0.062 ± 


0.140 


43, 


.652 ± 





.254 


SNLS- 


-04D3nh 





.340 


griz 


22, 


.137 ± 





.004 


1. 


.011 


± 0.006 


0.089 ± 


0.004 


41, 


.323 ± 


0, 


.012 


SNLS- 


-L)4D3nr 





.960 


griz 


24, 


.542 ± 





.075 


0. 


.922 


± 0.045 


0.070 ± 


0.110 


43, 


.622 ± 


0, 


.234 


SNLS- 


-04D3ny 





.810 


griz 


24, 


.272 ± 





.050 


1. 


.005 


± 0.084 


-0.065 ± 


; 0.152 


43, 


.691 ± 





.301 


SNLS- 


-04D3oe 





.756 


griz 


24, 


.069 ± 





.026 


0. 


.783 


± 0.028 


-0.259 ± 


: 0.033 


43, 


.453 ± 


0, 


.058 


SNLS- 


-04D4an 





.613 


griz 


24, 


.022 ± 


0. 


.023 


0. 


.823 


± 0.025 


0.064 ± 


0.025 


42, 


.961 ± 





.061 


SNLS- 


-04D4bk 





.840 


griz 


24, 


.314 ± 


0. 


.037 


1. 


.050 


±0.051 


0.142 ± 


0.098 


43, 


.475 ± 





.185 


SNLS- 


-04D4bq 





.550 


griz 


23, 


.362 ± 





.020 


0. 


.995 


± 0.029 


0.112± 


0.027 


42, 


.487 ± 





.056 


SNLS- 


-04D4dm 





.811 


griz 


24, 


.390 ± 





.044 


1. 


.000 


± 0.057 


-0.161 ± 


: 0.150 


43, 


.950 ± 


0, 


.264 


SNLS- 


-04D4dw 





.961 


griz 


24, 


.566 ± 





.093 


0. 


.962 


± 0.058 


-0.117 ± 


: 0.138 


44, 


.000 ± 





.290 



''Heliocentric redshift. 

•"See section lsm for description. 

'Computed with Hq = 70 km s 'Mpc^'. Uncertainty only accounts for photometric uncertainties. 
■•Not included in the final cosmological fits. 



